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Abstract 

We study Hermitian random matrix models with an external source matrix which has 
equispaced eigenvalues, and with an external field such that the limiting mean density of 
eigenvalues is supported on a single interval as the dimension tends to infinity. We obtain 
strong asymptotics for the multiple orthogonal polynomials associated to these models, and 
as a consequence for the average characteristic polynomials. One feature of the multiple 
orthogonal polynomials analyzed in this paper is that the number of orthogonality weights 
of the polynomials grows with the degree. Nevertheless we are able to characterize them 
in terms of a pair of 2 x 1 vector-valued Riemann-Hilbert problems, and to perform an 
asymptotic analysis of the Riemann-Hilbert problems. 

1 Introduction 

We consider random matrix ensembles under the influence of an external source matrix with 
equidistant eigenvalues. The ensembles consist of the space of n x n Hermitian matrices with a 
probability distribution of the form 

■^eicp{-nTr[V{M) - AM])dM, (1.1) 

where 

n 

clM = Y[d^M,jdQMij YldMjj. (1.2) 

i<j j=l 

The external field V{x) is a real analytic function which has sufficiently fast growth at infinity, 

lim Xi^ = +oo, (1.3) 

x-s>±oo |x| -|- 1 

and the external source matrix yl is a deterministic Hermitian matrix. Due to the unitary 
invar iance of the model, we assume, without loss of generality, 

^ = diag(ai,a2, . . . ,a„). (1.4) 

In our paper, we further assume the eigenvalues of A are equispaced on a certain interval, such 
that aj = a{j — l)/n + b where a and b are constants. By arguments of shifting and scaling, it 
suffices to consider the case 

j - I 

Oj = , where j = 1, 2, . . . , n, (1-5) 

n 
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and we work with the external source matrix A given by (|1.4p and p.5|) throughout this paper. 
The normahzation constant Z„ in (jl.ip depends on n and V. In the simplest example we have 

2 

V{x) = which gives the Gaussian Unitary Ensemble (GUE) in external source A. If we allow 
singularities of V, and let V{x) = (x — ■^logx)xx>0) we have the complex Wishart ensemble 
that has wide applications in statistics and wireless communication, see e.g. [8]. 

Random matrix ensembles with external source were introduced in [18[ [33] . and are inti- 
mately connected to multiple orthogonal polynomials [15]. If the external field is the classical 

2 

one V{x) = ^ or V{x) = {x — ^ logx)xa:>0) i-e-, the ensemble becomes the GUE with external 
source or the complex Wishart ensemble, more techniques are available for asymptotic analy- 
sis, and for a large class of external source matrices, including the equispaced one defined by 
()1.4p and (jl.Sp . the asymptotics can be obtained. See [22] for the complex Wishart ensemble. 
However, when the external field V{x) is general, the asymptotic analysis of the random matrix 
ensembles with external source has only had success for particular choices of external source 
matrices. Asymptotics for large n have been studied in [HI [16l [5l HZl [H |13l |6] in the case where 
the external source matrix A has two different eigenvalues a and —a with equal multiplicity, 
and in [9l [TOl \TT\ [T2] when A has a bounded, or slowly growing with n, number of non-zero 
eigenvalues. Large n asymptotics for general external source matrices have been considered in 
the physics literature, see e.g. [23], but rigorous asymptotic results have not been obtained to 
the best of our knowledge except for the two above-mentioned cases. We remark that the GUE 
with external source and the complex Wishart ensemble have other generalizations, the complex 
Wigner matrix model with external source and the complex sample covariance matrix model 
respectively. They have also been studied extensively, see e.g. [7]. 

Let us first recall some general properties about random matrix ensembles with external 
source. An ensemble of the form (jl.ip with eigenvalues of the external source matrix being 
ai, . . . , On induces a probability distribution on the eigenvalues Ai, . . . , A„ of the matrices given 
by [m [211 [27] 

n n 

— det(e""'^0^,=i A(A) J] ^-^^^'^^ [J '^^r (l-^) 

i=i i=i 

where Z'^ = const -Zn ■ A(a), and A(A) = rij<;j(Aj — Aj) and A(a) = Wi^j{aj — are Vander- 
monde determinants. A remarkable fact is that the average characteristic polynomial of such 
an ensemble (jl.ip satisfies orthogonality conditions: indeed, let 

n 

:= E„(det(zJ - M)) = K{\{{z - A,-)), (1.7) 

where E„ is the average with respect to (jl.ip . and is the average with respect to (jl.6p . then 
it was proved in [15| that p^'^ is characterized as the unique monic polynomial of degree n 
satisfying the orthogonality conditions 

/ pi")(x)e"'^^^e-"^(^)(ix = 0, for j = 1, . . . ,n. (1.8) 

These are the orthogonality conditions for the so-called type II multiple orthogonal polynomials 
with respect to n different orthogonality weights e^^J^e""^^^^ j = 1, . . . ,n. Specialized to our 
situation aj = for j = 1, . . . ,n, the joint probability distribution of the eigenvalues takes 
the form 

^ n n 

" i<j i<j i=l j=l 
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and the monic type II multiple orthogonal polynomials Pj {x), where j = 0, 1, . . . is the degree, 
are characterized by 

/p5"^(x)e^^e-"^(^)dx = 0, for A; = 0,...,j-1. (1.10) 

It is well-known that the point process (|1.6p is determinantal [33], and its two-point correla- 
tion kernel can be written in terms of multiple orthogonal polynomials. If aj = ^-^^ the kernel 
takes the form [15] 

n— 1 ^ 

Kn{x,y) = e-Wi^^^viv)) ^ -L.pf\x)Qf\yl (1.11) 

where p^J^\x) are the type II monic multiple orthogonal polynomials characterized by (|1.10p . 

and Q^j^\y) = (e^) are linear combinations of e^^ with k = 0, 1,2, ... ,j, subjected to the 
orthogonality conditions 

/x'=gf^(e^)e-"^(^)dx = 0, for A; = 0, . . . , j - 1, (1.12) 
Jr 

where gj"'* is a monic polynomial of degree j. Finally the constants h^"'^ are given by 

/iW = / p'-;'\x)qPie^)e-^^(^Ux. (1.13) 
Jr 

The orthogonality conditions (jl.lOp and ()1.12p for p^"^ and Qj^^ can also be written at once as 

/ pS"^(x)gi"^(e^>-"^(^)dx = 0, for j^keN. (1.14) 
Jr 

Note that the multiple weights g'^^e""^*^^^ constitute an AT system |3m Section 4.4], and hence 
p^"^ and Qj^^ are uniquely defined, and /i^"^ 7^ [3T] . 

Remark 1. As the counterpart of p^j^\x), Q^"'\x) is the j-th multiple orthogonal polynomial 

(n) 

of type I, up to the constant factor hj . Generally the type I multiple orthogonal polynomials 

are not polynomials, but in the present setting, Qj^\x) is a polynomial in e^. 

Remark 2. If the external field ^ is a quadratic polynomial, distributions of the form (jl.6p 
can also be realized in models consisting of n non-intersecting Brownian motions. In particular, 
(jl.9p is the joint probability distribution at an intermediate time of n non- intersecting Brownian 
motions starting at one point and ending at n equidistant points. Such a model has been studied 
in [29]. Different endpoint configurations have been studied e.g. in [21 [3]. 

In analogy to ()1.7p . qlC^ can also be interpreted as an average over the determinantal point 
process (|1.9p . We will prove the following result in Appendix lA.il 

Proposition 1. Let V be real analytic satisfying ()1.3p . We have the identities 

n 

ql^He^) = E„(det(e^^ - e^)) = KiU^e^ " e'O), (1-15) 

i=i 

where E„ denotes the expectation associated to (jl.ip with A given by ()1.4p ~ (|1.5p . and is the 
expectation associated to (jl.9p . 
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The main goal of this paper is to obtain asymptotics for the average characteristic polynomi- 

(n) 

als pn of the random matrix ensemble as n — )• oo. In addition we will also obtain asymptotics 

(n) 

for the dual polynomials g„ . A crucial role in the description of the asymptotic behavior of the 
polynomials will be played by an equilibrium measure. By ()1.9p . the joint probability density 
function of eigenvalues is maximal for the n-tuples (Ai, . . . , A„) for which 



log - Xj\-^ + ^log\e^' - e^^-^ + nY^ViXj) (1.16) 



i<j i<j j=l 



is minimal. As in [211 Section 6.1], one can then expect that the limiting mean distribution of 
the eigenvalues of the random matrices is given by the equilibrium measure /Uy which minimizes 
the energy functional 

Ivif^) ■=\jj log\t-s\-^d^l{t)d^x{s) + ]^ II log\e' -eT^dfi{t)df,{s) + ^ V{s)dfi{s), (1.17) 

among all Borel probability measures // supported on M. This is in analogy to the equilibrium 
measure corresponding to a matrix model of the form (jl.ip without external source, which is 
given as the unique minimizer of the energy 

log\t- s\-^dn{t)dn{s) + I V{s)d^j.{s). (1.18) 

Following the proof in [21] of existence and uniqueness of the minimizer of ()1.18p , we will show 
existence and uniqueness of the equilibrium measure minimizing (jl.l7p . 

Theorem 1. Let V he real analytic, satisfying the growth condition (|1.3p . Then there exists a 
unique measure fj, = fiy with compact support which minimizes the functional (jl.lTp among all 
probability measures on M. 

Remark 3. It should be noted that the growth condition (jl.Sp is stronger than the usual loga- 
rithmic growth needed to have a unique minimizer of the one- matrix logarithmic energy (jl.lSp . 
This is a consequence of the second term in (jl.lTp . 

The proof of this result will be given in Section [2] but it does not give any information about 
the measure itself. For that reason, in what follows, we will restrict ourselves to a class of 
external fields V for which the equilibrium measure behaves nicely and is supported on a single 
interval. 

We say that a real analytic external field V satisfying (|1.3p is one- cut regular if there exists 
an absolutely continuous measure d^y{x) = il)y{x)dx satisfying the properties 

(i) supp = [o-j for a < & € M, and / dfiv{x) = 1, 

(ii) ijjvix) > for X € (a, 6), 

(iii) lima;_s.a , ^nd \\m.x-^i)_ exist and are different from zero, 

(iv) for X G [a, 6], there exists a constant i depending on V such that 

I \og\t - x\-^dfj,v{t) + I log|e* - e'^l'^dfivit) + V{x) + i = 0, (1.19) 

(v) for X G M \ [a, 6], we have 

I log\t - x\~^dfivit) + I log|e* - e'^l-^dfivit) + V{x) + i > 0. (1.20) 
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Properties (iv) and (v) are variational conditions for fiy, and it follows from standard arguments 
that a measure satisfying (i) , (ii) , (iv) and (v) minimizes the energy functional ()1.17p . Under 
the condition that V is one-cut regular, we obtain large n asymptotics for p^\z) and qi^\e^) 
defined by ()1.14p . and state it in the following theorem. For the purpose of a subsequent paper, 
we give slightly more general asymptotics for ^^'^^(z) and (7jj^^(e^), where A; is a constant 
integer. 

Suppose the equilibrium measure fiy associated to V is supported on a single interval [a, b] 
and the density function is ipvi^)- In order to be able to formulate our results, let us define 
Co G M and ci G M"*" such that 



ci 



1 1 

7 + — 
4 ci 



Co 



log 



b + a 



4 + ci 



1 + i + 1 

4 ^ ci ^ 2 



(1.21) 
(1.22) 



Note that ci is well defined since as ci runs from to +00, the left-hand side of (|1.22|) increases 
monotonically from to -|-oo. Then we define the transformation 



J(s) = Jci,co{s) ■= cis + Co - log 



s + 



(1.23) 



for s G C \ [—2, 2]) where the logarithm corresponds to arguments between — vr and tt. For 
s < — ^, 3ci,co{s) has a maximum at Sa, and for s > ^, 3ci,co{s) has a minimum at S},, where 



1 1 

4 + ^' 



Sb 



1 1 

4 ci 



(1.24) 



The extrema Sa and Sfe are also characterized by identities a = Jci,co(sa) and b = Jci,co(sfe)- 
In Section [221 a region D C C is defined by Proposition [21 and it is shown there that J maps 

biholomorphically into S \ [a, 6], 



.1 il 

2' 21 



C \ D biholomorphically into C\[a,b], and maps D \ 
where 

S := {z G C I -TT < < vr}. (1.25) 

Let the functions Ii and I2 be inverse functions of J for these two branches respectively: Ii is 
the inverse map of Jco,ci from C \ [a, b] to C \ D, and I2 is the inverse map of Jco.ci from § \ [a, b] 

2' 2J' 



toZ)\[-i,i]: 



Ii(J(.)) 



s, 
s, 



for seC\D, 
iovseD\[-^,^]. 



(1.26) 
(1.27) 



Writing, for x G (a, 6), 



(x) :- 



lim Ii(x -|- ie) 

e-S>0+ 

lim Ii(x — ie) 

e-S>0+ 



lim l2{x — ie) 

e-S>0+ 

lim lofx -|- ie) 



(1.28) 
(1.29) 



we have that I+(x) lies in the upper half plane, I-(x) in the lower half plane, and their loci are 
the upper and lower boundaries of D (denoted as 71 and 72 in Proposition [2|) respectively. 
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Let the functions Gk{s) and Gkis) be defined as 

G„.):..f(i±i>(i^, <i^, (1.30) 



where the square root y s'^ — \ — ^ has its branch cut along the upper edge of D (71 defined 
in Proposition [2|) in Gk{s), along the lower edge of D (72 defined in Proposition [2]) in Gk{s), 
and \J s"^ — \ — ^ ~ s as s — )• c« in both cases. Further we define 

rfc(x) := 2|Gfc(I+(x))|, O^ix) := arg(Gfc(I+(x))), (1.31) 

ffc(x) := 2\Gk{I-{x))l Okix) := arg(Gfc(I_(x))), (1.32) 

for X in (a, 6). 

We also need to define the functions 

giz) := / log{z - x)ilJvix)dx, g{z) := / log(e'^ - e^)^/'y(x)dx, (1.33) 

./a ./a 

with the branch cut of the logarithms for z G (—00, x) and E (0, e^), and ■i/'v' is the equilibrium 
density. Let 

(l){z):=g(z)+g{z)-V(z)-i (1.34) 

for 2; G S \ (—00, 6), where £ is a constant to make </)(a) = (/)(6) = (see ()1.19p and ()3.8p ). Then 
we will see later on, see Section 14.41 that 

/,(z) := (-^^(z)^ (1.35) 

is a well defined analytic function in a certain neighborhood of b, with /b(6) = 0, fjj{b) > 0. 
Similarly, 

2 

fa{z) := (-^cA(z) ± ^TTz) ' (1.36) 

(where the sign is + in C"*" and — in C~,) is a well defined analytic function in a certain 
neighborhood of a, with fa{a) = 0, /^(a) < 0. 

Since bothpj"^(z) and qj^\e^) are real analytic functions, it suffices to give their asymptotics 
in the upper half plane and the real line. For the ease of the statement of the theorem, we divide 
the upper half plane into regions ^4^, Bg, Cs and where (5 is a small enough positive parameter, 
such that Cs and are semicircles with radius 5 and centered at a and b respectively, 
consists of complex numbers not in Cs or Ds, with real part between a and b and imaginary 
part between and |, and As = C+ \ {Bs U C5 U Ds). See Figured] for the shapes of the four 
regions. 

Theorem 2. Let V be one-cut regular. As n ^ 00, we have the following asymptotics of 
P'n+k^^) and g^'^^(e^), A; G uniform for z in regions As, Bs, Cs and Ds, if 8 is small enough. 

(a) In region As and on its boundary, 

plri(^) = (1 + 0(n-i))Gfe(Ii(z))e"s(^), (1.37) 
q'nlM) = (1 + 0(n-i))G,(l2(z))e"^(^), (1.38) 
where (jl.38p is valid for < vr. 
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a 



b 



Figure 1: The four regions in the complex upper half plane where asymptotics for the multiple 
orthogonal polynomials P^^+k^^) ^^'^ 1n+ki^^) ^^^^ given in different formulas. 



(b) In region Bs, 

Pnlk{z) = (1 + 0(n-i))Gfc(Ii(z))e"^(^) + (1 + 0(n-i))G,.(l2(z))e"(^(^)-i(^)+^), (1.39) 

. (1.40) 



= (1 + 0(n-^))Gfc(l2(z))e"g^^^ + (1 + 0{n-'))Gk{li{z))e 



in+k 



ng{z) 



If X £ (a, 6) is on the boundary of region Bs, then 



P^nlkix) = rfc(:E)e"/'°gl"-^l'^''^(^) cos Lt: j\iiv{t) + 6 k{x) \ + 0{n-^) 



fk{x)e 



n f\og\e''-ey\d^^v{y) 



cosfnvr / d^xv{t) + (>k{x)] + 0{n ^) 



(1.41) 
. (1.42) 



(c) In region Cs, let Ai denote the Airy function 'Jj. Then 



P%{z) = ^ ((1 + 0(n-i))Gfc(Ii(z)) - (1 + 0{n-^))iGk{Uz))) nl f^z) M{ni fa{z)) 
- ((1 + 0{n-^))Gk{li{z)) + (1 + 0{n'^))iGk{l2{z))) n-lfa~'{z) Ai'{nlfa{z)) 



qS-ki^n = V^ ((1 + Oin-'))GkiUz)) - (1 + 0{n-'))iGk{Iiiz))) J fliz) Ai(ni/„(z)) 
- ((1 + 0{n-'))Gk{l2{z)) + (1 + 0{n-'))iGkiIiiz))) n'T, ff^ {z) k\' {nl fa{z)) 

1 1 
where fa{z) has branch cut on {a,b), and fa{x) > for x < a. In particular, if z = 

a + f'a{a)~^n~'^^^t with t bounded, then 

p^(g{z)-s{z)-V{z)~e) in) ( . _ 



(_l)V27r - + 



1 1 



4 ci 



1 1 ^ fc-i 

7 + — + o 
4 ci 2 



c,'H-f'M)-'nl {k\{t) + 0{n'l)), 

(1.45) 



(-1)^^2^(1 + 1 



1 1 1 

4 ci 2 



-k 



M^+^^^){-f'^{a))^J (^Ai(t) + 0{n-i)^ . 



(1.46) 
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(d) In region Ds, 



((1 + 0{n-^))Gk{li{z)) - (1 + 0{n-^))iGk{l2{z))) f} {z) M{nl h{z)) 

((1 + 0{n-^))Gk{li{z)) + (1 + 0{n-^))iGk{l2{z))) n-lf^'^iz) Ai'(ni ^(z)) 

, (1.47) 



2ti\\ + - 
4 ci 



2%{el = ^ ((1 + 0(n-i))G'fc(l2(z)) - (1 + 0{n-^))iGk{Uz))) n\f^{z) Ai{J Mz)) 
\l + 0{n~'))Gk{l2{z)) + (1 + 0{n-'))iGk{Ii{z))) n-lf^^z) A\'{nlh{z)) 

If z = b + fl{b)^^n~'^/^t with t bounded, then 
et(i(^)-gW-V'W-^)pW^(,) = 



fc— 1 

i + l-i) c^-V^(?')^n^(Ai(t) + 0(n-t)), (1.49) 



Il(giz)-Siz)-V(z)~£) (n) ( z 



27r I 7 + — 

4 ci 



1 1 

4 ci 



1 \-*^ 



/^(6)4n6 Ai(i) + 0(n- 



(1.50) 



fej r/ie inner product of p^^]^^{z) and q^^Ji^ie^) defined in (jl.l3p /las ^/le asymptotics 

h^^l^ = 27rq+^e'=(^+^'')e«^(l + ©(n'^)). (1.51) 

In general it is not easy to determine whether an external field V is one-cut regular or not, 
or to find the support [a, b] of the equilibrium measure and the density function tpy However, 
if the external field is strongly convex, i.e. V"{x) is bounded from below by a positive constant 
for a; G M, then V is one-cut regular, and we can compute the support and density function of 
the equilibrium measure explicitly in terms of the functions I± defined before. 

Theorem 3. If V is a real analytic strongly convex function, then V is one- cut regular. More- 
over, the quantities cq and c\ that are related to the endpoints a, b of the support of the equilib- 
rium measure by ()1.2ip and ()1.22p are obtained by solving a pair of equations ()3.2p and (|3.3p 
expressed in V, and a,b are determined by ci and cq by (|3.4p . (|1.2ip and p.22p . The density 
function tpy is given by 



1 

2^ 



V"{u) log 



I+(7i)-I_(x) 



I+(n)-I+(x) 



du. 



(1.52) 



Remark 4. The conditions of Theorem [3] are sufficient but far from necessary to have one-cut 
regularity. See Example 2 in Appendix |B] for non-convex one-cut regular external fields. 

For the random matrix model without external source, it is well known that 
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(1) the empirical distribution of the eigenvalues of the random matrix, 



(2) the normalized counting measure of the n-Fekete set, 

(3) the normalized counting measure of the zeros of the orthogonal polynomial (which is the 
average characteristic polynomial of the random matrix), 

all converge to the equilibrium measure as the dimension n — )• oo. The counterpart of |(3)| in 
our equispaced external source model, in case that the external field V is one-cut regular, is a 
direct consequence of Theorem ^[h)\ 

Corollary 1. Let V be one- cut regular, and and q>a'^ be defined by p.lOp and p.l2p 

respectively. Suppose real numbers Zj and Zj are zeros of p^\z) and q^\e^) respectively, and 
l^n = ^ Yl]=i and pLn = ^ 5]^j=i respectively. Then as n ^ oo, fin o,nd fin converge weakly 
to /iy. 



Counterparts of (1), (2) and |(3)| can also be proved by mimicking the arguments in 



Sections 6.3 and 6.4]. Although we are not going to pursue this approach, we remark that all 



the counterparts of (l) -(3) should not rely on the assumption of one-cut regularity. 



Outline 

In Section [21 we prove the uniqueness and existence of the equilibrium measure, as stated in 
Theorem [TJ In Section [3l we explain in detail how one can construct the equilibrium measure 
and its density in the case of a strongly convex external field V , by solving a scalar Riemann- 
Hilbert (RH) problem and by using the transformation J. This also leads to the proof of 
Theorem [3l In Section U we characterize the polynomials p^^+k ™ terms of a 1 x 2 RH problem. 
This RH problem is different from the known [n + k + \)y.{n + k + \) RH problem characterizing 
the multiple orthogonal polynomials p^^+k ES] and from the classical RH problem for orthogonal 
polynomials [23]. Since re is a large parameter in our settings, the 1x2 RH problem is much 
more convenient for asymptotic analysis. As a drawback, our RH problem is non-standard in 
the sense that the entries of the solution live in different domains. This requires a modification 
of the Deift/Zhou steepest descent method to analyze the RH problem asymptotically. Again 
the transformation J plays a crucial role here. It allows us to transform the 1x2 RH problem 
to a scalar shifted RH problem, and we then obtain small norm estimates for the solution to 
this RH problem. In Section [5l we formulate a similar RH problem and perform a similar 
asymptotic analysis for the polynomials q^^j^ • In Section [6l we use the results obtained from 
the RH analysis to prove Theorem [2] and Corollary [TJ In Appendix [XJ we prove Proposition [1] 
and several technical lemmas used in this paper. In Appendix [Bj we give explicit formulas for 
the equilibrium measure for quadratic and quartic V as examples. In Appendix[Owe derive the 

(n) 

asymptotics for the polynomials p)i for quadratic V using an integral representation and the 
classical steepest descent method. In this derivation we show that the transformation J also 
arises in a more direct way in the equispaced external source model. 

The main novel feature of this paper is the successful asymptotic analysis of the non-standard 
RH problem which characterizes the multiple orthogonal polynomials. Although the resulting 
large re asymptotics for the polynomials resemble those for usual orthogonal polynomials rel- 
evant in the one-matrix model without external source, the RH method used to obtain those 
asymptotics had to be modified in a nontrivial way. We feel that the modification of the RH 
method, with in particular the use of the transformation J, is the main contribution of the 
present paper. We believe it is the first time that a RH analysis has been carried through for 
multiple orthogonal polynomials with a growing number of orthogonality weights. 
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2 Proof of Theorem [T] 

Following [211 Section 6.2] (see also [2E])) one can prove the existence of a unique Borel proba- 
bility measure minimizing the energy /y(^) given in p.lZp . which can conveniently be written 
as 

Iv{l^)= [[ k^it,s)df,{t)dfi{s), (2.1) 



with 

k'^it, s) = ^ log \t-s\-' + ^ log |e* - eT' + ^Vit) + ^V{s). (2.2) 



From the inequality |v — u| < Vl + f ^Vl + u'^ for u G M, we obtain 

Uog\t-s\-' + ^log\e'-eT'>-\iog{l + t^)-^log{l + s^)-^log{l + e^')-^log{l + e^n- 

(2.3) 

If V satisfies the growth condition (jl.3|) . it easily follows that there exists a constant cy such 
that {t, s) > Cy for all s, t G R. Thus /y (/x) > cy for any probability measure //, which implies 
that Ey = inf{Iy(/i)} > cy, where the infimum is taken over all probability measures on M. 
This is the crucial estimate for proving the existence of a unique equilibrium measure. The 
existence follows, exactly as in [21| Section 6.2], from the construction of a vaguely convergent 
tight sequence fin of measures with limit such that /y(/i) = Ey, as well as the fact that any 
minimizer must have compact support. 

The uniqueness is slightly more complicated, and we need the following lemma for it: 

Lemma 1. Let fi be a finite signed measure on M such that J dfi = and with compact support. 
Then 

\og \x - y\^^dfi{x)dfi{y) > 0, (2.4) 

log |e^ - ey\-Un{x)dfx{y) > 0. (2.5) 

The first inequality was showed in |2H Lemma 6.41], and the second part can be proved by 
replacing x i— )• e^' and y i— t- in the proof. 

Now assume that we have two measures //y and fl such that Iy(/iy) = /y(/i) = Ey- Then, 
for nt = fly + t{jl — fj,y) and t G [0, 1], we have 

1 1 ~ f 

Ivil^t) = ^I{fJ-v,fJ'v) + ^Hf^v^^J-v) + / V{x)dfiv(x) 



+ t yl{fiv,fl - fiy) + I{fiy,il - My) + J V{x)d{jl - Hv){x)) 

+ - MF, A - ^J'V) + ^^(/i - /iy, /i - ^y) j , (2.6) 

where 

/(/X, v) = jj log \x - y\-^dii{x)du{y), (2.7) 

/(m, v) = 11 log - ey\-Ufi{x)du{y). (2.8) 



10 



The above lemma ensures that lyipt) is a convex function of t. But since /it is a probabihty 
measure, we have IvifJ-t) > IvifJ-o) = IvifJ-i) = Ey, and hence /y(/it) = Ey for aU t G [0, 1]. In 
particular this implies 

- fiv,fJ'- fJ'v) + ^HP- - f^v,fj'- f^v) = 0, (2.9) 

and using a similar argument as in [21J, this implies that //y = ft, which yields the uniqueness 
of the equilibrium measure. 



3 Construction of the equilibrium measure 

In this section we assume the external field y is a convex real analytic function and V"{x) is 
bounded below by a positive constant for all x G R. We are going to show that V is one-cut 
regular, by an explicit construction of its equilibrium measure. The strategy of our construction 
is as follows. First in Section [3.11 we give the support of the equilibrium measure [a,b] without 
proof. Then in Section 13.21 we compute the density of the equilibrium measure, based on the 
information of the support. The density function is expressed in terms of the so-called g- 
functions g{z), g(z) and their derivatives, which are characterized by a RH problem. At last in 
Section 13.31 we verify that the measure with the support and the density obtained in the first 
two steps satisfy the criteria of one-cut regularity, and conclude that it is the unique equilibrium 
measure that we want to construct. 

Remark 5. In what follows, it may seem that the values of the endpoints a and b appear out of 
the blue, but if the external field V{x) is quadratic, the endpoints (as well as g{x) and g(x)) can 
be computed by a classical steepest-descent method. This computation is shown in Appendix 
Oas our inspiration. 

Remark 6. If an external field is non-convex but we know a priori that it is one-cut regular with 
support [a,b], then the method in Section [3^2] can still be applied and allows us to obtain the 
expression of the density function of the equilibrium measure. 



3.1 The support of the equilibrium measure 

Let 

Jxi,xo be defined as before by 

Jxi,xo{s) = XiS + xo -log ^, (3.1) 

and let 7 = Jxi,xo(i'^^^])^ depending on xi,xq, be the boundary of the region D defined in the 
Introduction, consisting of the curves 71 and 72, encircling the interval [— ^, |] in the counter- 
clockwise direction, see also Proposition [2] below. Since Jxi,xo{s) £ [a, &] for s G 7, V' {Jci^coi^)) 
is well defined for s in a neighborhood of the curve 7, if V is real analytic. 

Lemma 2. Given any strongly convex real analytic function V, i.e. such that V"{x) > c > 
for all X GM, the system of equations with unknowns xq and xi 

= 1^. iv'{3x,,x,{s))ds, (3.2) 
has a solution xq = cq G M and xi = ci G M"*". 
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We will prove Lemma[2]in Appendix A. Based on this lemma, we construct the support, and 
furthermore the density function, of the equilibrium measure. We do not prove the uniqueness of 
the solution of equations i\'6.2\i and {\'6.'6^ , for this uniqueness is a consequence of the uniqueness 
of the equilibrium measure by Theorem [H as from different solutions we construct different 
equilibrium measures. 

Here and later we take the value of the parameters cq and ci as the pair of solutions of ()3.2p 
and p.Sp . Then we claim that a and b, the two edges of the support of the equilibrium measure, 
are given by 

a = Jci,co('Sa), 6 = Jci,co(Sfe), (3.4) 

where Sa = — y^^|~T~^) Sb = \ + Then it is easy to verify that equations (|1.2ip - ()1.22p are 
satisfied. 



3.2 The g-functions and the density function of the equihbrium measure 



Under the assumption that the external field V is one-cut regular, with equilibrium measure 
diiv{x) = ijjv{x)dx supported on [a,b] as we claimed in ()3.4p . we construct two functions 
g{z) = /log(z — x)d^y{x) and g,{z) = J"log(e^ — e^)d^v{x) as in ()1.33p . To describe the 
domain of the function e^^^\ we introduce the notation of the cylinder S*^ which is formed by 
identifying the two edges of the strip §. If a function f{z) is defined for z E §, the limits 
/(x ± vrz) = \\m.z^x±TTi,z& f{z) exist point- wise, and furthermore f{x + ni) = f{x — iri), we say 
/ is defined on 'B'^. The properties (i) - (v) in the Introduction satisfied by nv are then translated 
into properties satisfied by g and g as follows. 

(i) For X G (— oo, a), 

g+(x) = g_(a;) + 27ri, g+(x) = g_ (x) + 27ri, (3.5) 

and then e^{z) is analytic in C\ [a, b] and e^*^^-* is analytic on the cylinder with slit [a, b]; 
es(^) ~ z as z — > oo, e^^^^ ~ as ^z — > +oo and e^*^^) = 0(1) as ^z -oo. 



(ii) for X £ (a, b), we have 

Tpv{x) = 



--^(g'+(x)-g'_(x)) = --^(gV 

ZTTl ZTTl 



[x) - gLix)) > 0, 



(3.6) 



(iii) as z — )• a or z — )• 6, the limits of g{z), g{z), g'{z) and g,'{z) exist, and as x — )• a+ or x — )• 6_ 
for X € (a, b), 



y - g'-{x) gV(x) 
iim . , hm 



X — a 



X — a 



lim 



sV(x) - g^_(x) 
y/b — X 



lim 



all exist and are all different from zero, 

(iv) for X £ [a, b], there exists a constant i such that 

g±(x) + g=p(x) - V{x) 

(v) for X G M \ [a, 6], we have 

g±(x) + g^(x) - y(x) 

Let us consider the derivatives 

G{x) := g'(x) and G{x) := g'(x). 



0, 



< 0, 



gV(x) - g:.(x) 

y/b — X 
(3.7) 



(3.8) 
(3.9) 
(3.10) 



The properties (i) , (iii) and (iv) for g(x) and g(x) then imply that G and G need to satisfy the 
following RH problem: 
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RH problem for G and G 

(a) G is analytic in C \ [a, 6], G is analytic in S'^ \ [a, 

(b) for X G [a, 6], we have 

G±(x) + G'^(x) = 0, (3.11) 

(c) we have the asymptotic conditions that G±{x) and G±{x) are bounded for all x G [a, 6], 
and 

G(z) = ^ + C'(z"2), asz^oo, (3.12) 

G{z) = 1 + 0(6^^), as 3f?z ^ +00, (3.13) 

G{z) = 0(1), as -oo, (3.14) 

The main technical difficulty in solving the RH problem for G and G lies in the fact that 
the two functions live on different domains: G is defined in the complex plane with slit [a,b], 
and G is defined in the cylinder S"^ with slit [a,b]. In order to resolve this problem, we will use 
the transformation (jl.23p J(s) that maps C \ [— ^, ^] to both C and §. Recall that Sa and 
are the two critical points of J(s) given by ()1.24p and that they satisfy the identity ()3.4p . The 
following property will be used in the construction of G and G. 

Proposition 2. There are an arc 71 from Sa to Sb in the upper half plane, and an arc 72 from 
Sa to Sf) in the lower half plane, such that 

(a) J(7i) = J(72) = [o-,b], and the mapping is homeomorphic on these two curves. 

(h) Denote the region enclosed by 71 and 72 by D. Then J{C\D) = C\[a,b], and the mapping 
is univalent. 

(c) 3{D \ [— ^, ^]) = S \ [a,b], the mapping is univalent, and the upper and lower sides of 
(— ^, ^) are mapped to M — vri and M + vri respectively. 

The mapping J outside and inside D is illustrated in Figures [2] and [3l The proof of Propo- 
sition [2] is given in Appendix lA. 21 In Figured] we give two examples of 71 and 72 by numerical 
simulation. 




b 



Figure 2: Mapping J outside of D. 



Let us now define the function M{s) by 

M( forsGC\A .3^5. 

\g(J(.)) for.Gl)\[-l,i], 

so that M is analytic in C \ (71 U 72 U [— ^ , ^] ) . Note that the domain of G can be extended from 
§ to S"^, so that M{s) can be analytically continued to (— ^, ^) accordingly. The RH conditions 
for G, G are now transformed to the following conditions for M. 
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Figure 3: Mapping J inside of D. 







0.6- 

/ 0.4 - 




/ 0.2 - 




[ -1 -0.5 
\ 0.2 


^^^^^^^^^^ 


\ -0.4 - 
-0.6 - 




^V^^ -0.8 - 







Figure 4: The shapes of 71 and 72 when c\ = \ (left) and c\ = 10 (right), where 71 is the upper 
boundary of the region and 72 the lower boundary. 



RH problem for M 

(a) M is analytic in C \ (71 U 72 U {-\, i}), 

(b) M satisfies the jump condition 

M+(s) + M_(s) = V' (J(s)) , for s G 71 U 72, 

(c) M±(s) is bounded on 71 and 72, and M has the asymptotics 



1 

2' 



M(s) = — + 0(s"^), as s ^ 00, 

M(s) = l + 0(s- ^), ass 

M(s) = 0(1), as s ^ - 
It is straightforward to solve this scalar RH problem. We write 

[/(.) = y'(j(.)), 



(3.16) 

(3.17) 
(3.18) 
(3.19) 

(3.20) 



and note that U is analytic in a neighborhood of 71 U 72, since V is real analytic. Then it is 
readily verified that the unique solution M to the above RH problem for M is given by 



M{s) 



d^, for s G C \ D, 



for s e D, 



(3.21) 
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where 7 is the closed curve which is the union of 71 and 72 and has counterclockwise orientation. 
In particular, ()3.17p and (|3.18p follow from the system of equations ()3.2p and (|3.3p in Lemma E] 
satisfied by cq , ci . 

Now we can give an expression for g{z), g{z) and the density function ipvi^) of the equi- 
librium measure, under the assumption that the support of the equilibrium measure is known. 
Recall that Ii is the inverse map of J from C \ [a, 6] to C \ -D, I2 is the inverse map of J from 
C \ [a, b] to D, and their boundary values define I±(x), see (|1.26p - (jl.29p . We have I+(x) € 71, 
I_(x) € 72, and = To obtain a formula for the density ^y(x) of the equilibrium 

measure, note that it follows from (|3.6p and the identities G = g' ,G = g' that 



^v{x) = - G.{x)) = -^{G+{x) - G.{x)), for x G [a,b] 

From and (IXT^ . we obtain 



(3.22) 



^v{x) = -^(M+(I+(x)) - M_(I_(x))) = -^(M+(I_(z)) - M_(I+(x))), for x G [a, h]. 

(3.23) 

Applying the first identity in (|3.23p and the formula (|3.2ip for M(s), we let z = a; + ie, e > 0, 
approach x from above and have 



hm 



u{0 



1 



1 



.™47r2y, ^""^ l^I+(u)-Ii(z) I+(u)-Ii(z) I-(n)-Ii(z) I_(n)-Ii(z) 



lim — 77 



lV(n) 



lV(n) 



I+(n)-Ii(z) I+(n)-Ii(z) 
lim^ /"V'(7x)3f?-^log I 
1 



I+(^) -Ii(^) 
-Ii(^) 



du 



lim 



V>)Kiog[H^^l:iiMun 



^I+(^x)-Ii(z)^ 

iM-l-ix) 



U{u)-U{x) 



du. 



du 



(3.24) 



3.3 Proof of Theorem H 

We showed so far that the equilibrium measure associated to the external field V has the density 
function tpy as we have constructed in Section [3T2l as long as it is supported on the single interval 
[a, 6] that is given by ()3.4p . However, we have not proved that [a, 5] is the correct support yet. 
We will show that the measure with support [a, 6] and density function ijjyix) satisfies the 
properties (i) - (v) stated in the Introduction for one-cut regular equilibrium measures, which 



implies that the constructed measure is indeed the true equilibrium measure. Note that these 



properties are equivalent to properties (i) - (v) in Section 13.2 



Prom the construction of ipv{x), it is normalized, i.e., ipv{x)dx = 1. This follows from 
the asymptotics of G and G, given in (|3.12p and (|3.13p . and the definitions ()1.33p of g and g, 
the antiderivatives of G and G. 
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For X G (a, 6), it is geometrically obvious that |I+(u) — > |I+(n) — I+(x)|, and then 

Klog((I+(ti) — — I+(x))) > for all u S {a,b). Substituting this inequality into 

(j3.24|) and noting that V" is positive, we have that tpvi^) > for all x £ (a, 6). Similarly we 
have tpvix) — )• for X — )• a+ and x — )• 6_. 



The identity (jl.lOp that gives condition (iv) in the Introduction, or equivalently the identity 



()3.8p that gives condition (iv) in Section \3.2\ is obvious from the construction of ipy Thus we 
only need to prove the remaining two properties for the equilibrium measure hold, i.e., il^v{x) 
vanishes like a square root as x — )• a+ or x — )• 6_, and G+(x) + G_(x) — V{x) < £ for x < a or 
X > b. 

Let the function H be defined by 

H{z) = (g{z) + G{z) - V'{z)) ^ , (3.25) 



It is well defined where G, G, V are defined, and it can only be discontinuous on [a, b]. However, 
by (IMI) and 

H+{x) = (G+ - G-f = -ATT^-^vixf = (G_ - G+f = H^{x). (3.26) 

Hence H{z) can be defined on (a, b) so that o, b become isolated singularities. If we express 
G{z) and G{z) in terms of M(s) and then by the contour integral as in (j3.15p and ()3.2ip . we 
find that G{z) and G{z) grows at most logarithmically at a and b. Thus a and b are removable 
singularities of H{z), and H{z) can be defined analytically in S where V is defined, i.e., an open 
region containing the real line. Furthermore, by (|3.26p and the fact that il^vix) — )• as x — t- a+ 
or X — )• 6_, we have that H{a) = H{b) = 0. 

To show that ipvi-i^) vanishes like a square root at a and b, by (j3.26p it suffices to show 
that a,b are simple zeros of H{z). We consider b first. From ()3.26p and ()3.25p . we see that 
H[x) changes sign as the real variable x increases around 6, so if b is not a simple zero, it has 
multiplicity at least 3, and then -^^JH{x), which is well defined for x G (6, oo), would tend to 
as X — )• 6+. But we have for all x > 6 

^G(x) + G{x) - V'{x)] = g"(x) + g'{x) - V"{x) 



dx 



b 



1 



V'H^) ( + J^^Zr^2 ) ds - V"{x) < -V"{x). 

(3.27) 



Since V"{x) is bounded below by a positive constant, ■^y^H{x) cannot approach 0. Thus b is 
a simple zero of H{z). Similarly a is a simple zero. 

To show that G+(x) + G'_(x) — V{x) < ^ for x > 6, we need only that G+(x) + G'_(x) — V{x) 
is decreasing, since at x = 6 the identity G+(x) + G_(x) — V{x) = i holds. The decreasing 
property is given by the negative derivative shown in ()3.27p . Similarly we can show that 
G+(x) + G_(x) - V{x) < £ for X < a. 

Now we have proved that the measure '^^yix) on [a, 6] satisfies all the properties for one- 
cut regular equilibrium measures, so it is the unique equilibrium measure associated to V . 
Combining the results we have obtained in this section, we prove Theorem [3l 



4 Asymptotic analysis for the type II multiple orthogonal poly- 
nomials 

In this section, we write {x), the monic multiple orthogonal polynomials of type II satisfying 
orthogonality relations (jl.lOp . as Pj{x) if there is no confusion. 
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4.1 RH problem characterizing the polynomials 

Recall that the j-th degree monic polynomial Pj{x) = p^"'\x) is characterized by the orthogo- 
nality (jl.lOp . Consider the following modified Cauchy transform of pj: 

which is well-defined for z G S \ M. Since e~"^^^-* is real analytic and vanishes rapidly as 
X — )• ±00, for any polynomial p(x), we have the following asymptotic expansion for Cp{z) as 
z € S and 5Rz — )• +cxd: 

Cp{z) = ^ [ P(^) e-"^(")dx 

. M (4.2) 

fc=0 ^-^^ ' 
for any M S N, uniformly in Thus due to the orthogonality, 

_^(n) 

C'p.(z) = —^e-^^'+i)" + 0(e-(-''+2)^), (4.3) 

ZTTZ 

where /z^"^ is given by ()1.13p . For 2; G M, a residue argument shows that 

(Cp,)+(^) - (C'pfe)-(x) = p,(x)e-"^We-^ (4.4) 
Hence we conclude that if we consider Pj{x) and Cpj{x) together and write them in vector form 

Y{z) = Y(^^^\z) := {pj{z),Cp,{z)), (4.5) 

they satisfy the conditions 
RH problem for Y 

(a) Y = (Yi, Y2), where Yi is an analytic function defined on C, and Y2 is an analytic function 

on \ M, 

(b) Y has continuous boundary values Y± when approaching the real line from above and 
below, and we have 

y+(x) = y_(x)f ^ J, forxeM, (4.6) 
(cl) as 2; — )• 00, Yi behaves as Yi{z) = z^ + ©(z^"^), 

(c2) as 00 (i.e., ^z +00), Y2 behaves as Y2{z) = C'(e~(-'+^)^); as e"^ (i.e., 
^z — )■ —00), Y2{z) remains bounded. 

Conversely, the RH problem for Y has a unique solution given by (j4.5p . We give a proof 
of the uniqueness of the RH problem for Y based on the uniqueness of the multiple orthogonal 
polynomials pj. 

Theorem 4. The solution to the RH problem for Y above has a unique solution, given by 
Yi{z) = pj(z) andY2 = Cpj{z), where pj{z) is the monic multiple orthogonal polynomial of type 
II defined by (jl.lOp . and Cpj{z) is given in (|4.ip . 
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Proof. First, (j4.6|) in the jump condition |(b)| implies that Yi is an entire function, and condition 



(cl) imphes that Yi grows like as z ^ oo. So Yi =: p is a monic polynomial of degree j. 



Now we show that if y = {Yi,Y2) satisfies all the conditions (a) - (c2) of the RH problem, 
then Y2 is given in terms of Yi = p by 



1 f p{s 

^2\Z) = 7 

By condition |(b)[ Y2 satisfies 



y^i^) = / ^e--y^^^Us. (4.7) 



Y2,+{x) - Y2,-.{x) = p(x)e-"^(")-^ (4.8) 

Consider the function 

U{u) = Y2{\ogu) [ Zl(^e-"^Wds, (4.9) 

where we take the principal branch of the logarithm with branch cut on R^. Obviously U{u) is 
analytic for u £ C\M. By the jump condition Y2 on the real line given by (|4.6p and the property 
that Y2{x + Tri) = Y2{x — ni), we verify that U+{u) = U-{u) for u G (0,oo) or n G (— oo,0), so 
that U is an analytic function for n G C \ {0}. Note that since p is a polynomial and e""^*^**) 
vanishes rapidly as s — )• itcx), we have 



I ^(£Le-nV(s)^<, ^ 0(^-1) asTx^oo. (4.11) 

Jr - -u 



1 j' p{s 



From ()4.10p we find that is a removable singularity of U (u) and then U (u) is an entire function. 
Then from (|4.1ip we have U{u) = by Liouville's theorem. Therefore (|4.7p is proved. 

At last we apply the expansion ()4.2p for M = j — 1 to Y2 given in ()4.7p . We see that the 
asymptotic condition Y2 = 0{e~^^~^^^^) implies that 



/ p(2;)e'=^e-"^(^)dx = 0, k = 0,...,j-l. 
Jr 



(4.12) 



Comparing this with (|1.10p . we see that p = Yi is indeed the monic multiple orthogonal poly- 
nomial pj. □ 

Below we take j = n+k where k a constant integer, and our goal is to obtain the asymptotics 
for Y = y("-+'='") as n — ;> 00. 

4.2 First transformation Y t-^ T 

Recah g{z) and g{z) defined in (fOHD on C \ (-00,6] and S \ (-00,6]. Denote Y = y("+'=.") 
and define T as follows: 

r(z) := e-f y(.) ("'"J^^^ .nL)^"^"' (4-13) 

where £ is the constant appearing in (jl.lOp and (j3.8p . and 173 = (0^*1)- Then T satisfies a RH 
problem with the same domain of analyticity as Y, but with a different asymptotic behavior 
and a different jump relation. 
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RH problem for T 

(a) T = (Ti, T2), where Ti is analytic in C \ M, and T2 is analytic in S'^ \ M, 

(b) T satisfies the jump relation 

T+{x) = T^{x)Jt{x), for x G M, (4.14) 

with 

/^■n{g-{x)-s+{x)) ^n{g-{x)+g+{x)~V{z)~e)-x\ 

Mx)=(^ en(g+(x)-g_W) J> (4-15) 

(cl) as z ^ 00, Ti behaves as Ti{z) = + C'(z~('''+^)), 

(c2) as ^ 00, T2 behaves as r2(z) = ©(e-^'^+i)^), and as 0, Ta behaves as T2 = 0(1). 
4.3 Second transformation T S 

For X G M\ [a, b], it follows from the analyticity of and (|3.9p that the jump matrix Jt{x) tends 
to the identity matrix exponentially fast in the limit n — >• cx). For x G (a, 5), we decompose the 
jump matrix into 

/I 0\ / ^n{^_{x)+g+{x)-V{x)-tj-x 

Jt{x) = {_n(l,-{x)+x [^n{-g+{x)-g^{x)+V{x)+e.)+x q 



where the function ^(z) = g(2;) + g(z) — V^(z) — ^ is defined as in (|1.34p . The function (t){x) has 
discontinuity on (—00,6]), and by ()3.5p and (jS.Sp it satisfies 

<t)+{x) = (j)-{x) + 47ri for x < a, (4-17) 

<j3+{x) = - <j3-{x) for x £ {a,b). (4.18) 




Figure 5: The lens Sg. 

Then we "open the lens", where the lens is a contour consisting of the real axis and 
two arcs from a to b. We assume that one of the two arcs lies in the upper half plane and 
denote it by Si, the other lies in the lower half plane and denote it by S2, see Figure [5l We do 
not fix the shape of S5 at this stage, but only require that S5 is in S and V is analytic in a 
simply-connected region containing S5. 



Define 



Siz) 



Tiz) 
T{z) 

T{z) 



1 

-n(l){z)+z 
1 

p-n4i{z)+z 



outside of the lens, 

in the lower part of the lens, 

in the upper part of the lens. 



(4.19) 
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From the definition of S, we see that S is discontinuous on the upper and lower arcs with jump 
matrix ( ^) . On (a, b), it follows from ()3.8p and ()4.16p that the jump matrix for S takes 
the form ( ). Summarizing, we have the following RH problem for S. 

RH problem for S 

(a) S = (81,82), where Si is analytic in C \ S5, and 82 is analytic in S'^ \ S^, and S5 = 
M U Si U S2 is the contour depicted in Figure [5l 



(b) we have 

8+{z) = 8-{z)Js(z), for z G S5, 
where (note that e'^^^-' is well defined for z G (—00, a) by (j4.17p ) 



(4.20) 



1 0' 

e-'\ 

-e^ y ' 

'l ^n<f>{z)-z\ 

1 J 



for z G Si U S2, 
for z G (a, b), 
for z G M \ [a,b]. 



(4.21) 



(cl) as z ^ 00, Si(z) = z'' + 0{z^'^), 

(c2) as — )• c«, 5*2 behaves as 82{z) = 0{e'^''~^^^^), and as — 0, 5*2 behaves as 5'2(z) = 
0(1). 

By ()3.8p . we have, for x G (a, b), 

</.'±(x) = g^(x) + g'±(x) - y'(x) = g'±(x) - g^(x) = ^27r#i/(x). (4.22) 

Since ipvi^) > for all x G {a,b), it follows from the Cauchy-Riemann conditions that 

sR0(z) = (g(x) + g{x) - V{x) -£)>0 (4.23) 

on both the upper arc and the lower arc, if these arcs are chosen sufficiently close to (a, 6). As 
a consequence, the jump matrices for 8 on the lenses tend to the identity matrix as n — t- 00. 
Uniform convergence breaks down when x approaches the endpoints a and b. We need to 
construct local parametrices near those points. 

4.4 Construction of local parametrices near a and b 

Define 

y, := y,{0 = J Ai(a;^C), for j = 0, 1, 2, (4.24) 



where w = e 3 and Ai is the Airy function. 
Let 



r := e"~M+ Ue~]R+ u: 



(4.25) 
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Figure 6: The contour T and the jump matrices for A. 

be the contour consisting of four rays oriented each from the left to the right shown in Figure 
[6l and define the 2x2 matrix- valued function ^ in C \ F as 



27re 4 



2ne 4 



2/0 -?/2 \ 

,yo -y'2) 

-2/1 -2/2 



for < arg C < ^ 



-2/'i 



for \ < arg C < vr, 



27re 4 



zi / -2/2 2/1 



(4.26) 



27re 4 



-2/2 2/'i, 
i / 2/0 2/n 

.2/^, y'J 



for — TT < arg ( < 
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' 3 ' 



for < argC < 0. 



Using the identity 2/0 + 2/1 + 2/2 = 0, the fact that the Airy function is an entire function, and the 
asymptotics as C —^ 00 of the Airy function, one verifies that A satisfies the following model RH 
problem. This RH problem (and equivalent forms of it) appeared many times in the literature 
and is often referred to as "the Airy RH problem" , see for example \20\ [2T| . 



RH problem for A 

(a) ^ is a 2 X 2 matrix- valued function analytic in C \ F. 

(b) A satisfies the following jump relations on F, 

'1 1' 



A+(C) = ^-(C) 



1 



A+{C) = A^{C)(\ J 



^+(C) = A_(C) 



1 
-1 



for arg C = 0, 
for arg Q = — or arg Q 



for arg = tt. 



-27r 



(c) A has the following behavior at infinity. 



1 

V2 



r^'^^(_\ })e-x-3(J + o(^^/2))e-K''''^^ asC^oo, 



uniformly for ( G C \ F. 



Using the regularity condition which says that lima, exists and is positive, and the 
formulas of g{z) and g{z), and noting in addition that (j){b) = 0, we obtain the following local 
behavior for (j) near b, 

(piz) = -c{z - bf/"^ + 0{\z - 6|^/2), as z ^ 6, where c> 0. (4.31) 



(4.27) 
(4.28) 
(4.29) 

(4.30) 
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Then in a neighborhood Ub of b, there is a unique analytic function fj, satisfying = 0, 

flib) > and 

lMzf/' = -lHz). (4.32) 

Now we choose the lens of in such a way that fb{z) maps the jump contour f/fcn for S 
on the jump contour T for A, and we define the 2x2 matrix-valued function P^^\z) on Ub \ 
as 

P^^\z) :=^(n2/3/^(z))e-3MW-^)'^3_ (4^33) 
Using the jump relations (|4.27p - (|4.29p for A and ()4.17|) and (|4.18|) for ^(z), one verifies that 

Pf (z) = P^_^\z)Js{z), for z G C/5 n (4.34) 

where J5 is given in (j4.2ip . Since the determinant of A is identically equal to 1, ^ is invertible, 
and so is p(^)(z) ior z £ UbHT^s- By and (220]), we have 

S+{z)p'l'\z)-^ = S-{z)P^^\z)-\ for zeUbH^s- (4.35) 

Similarly, near a, 

(piz) = -c{a- zf^'^ + 0{\a- zf/'^)±2Tri, as z ^ a, c > 0, (4.36) 

where the sign of ib27ri depends on whether z is in the upper or lower half plane. In a neigh- 
borhood Ua of a, there is a unique analytic function fa satisfying /^(a) < and 

lfa{zf^'' = -\Hz)±7Ti. (4.37) 

Again we can choose the lens S5 in such a way that fa{z) maps the jump contour Ua H S5 for S 
on the jump contour F for A. Then define the 2x2 matrix- valued function P^^-^z) on Ua \ S5 
as 

P(")(z) :=^(n2/3/„(z))e-5("'/'W"^)'^3, (4_38) 
Similarly to (|4.34p and ()4.35p . we have 

S+{z)Pi^\z)-^ = S.{z)p[''\z)-\ for z G n Sg. (4.39) 

Remark 7. Usually, a local parametrix serves as a local approximation to the solution of the RH 
problem. Since S is vector-valued and our local parametrices P^""^ and P^''^ are 2 x 2-valued, 
this is not quite true in our situation, but it will turn out later that large n asymptotics for 
S near a and b can be expressed in terms of P^^^ and P^*^), and thus in terms of the Airy 
function. Later in Section 14.61 we will build a vector- valued "global parametrix" p(°°) , which 
approximates S away from the endpoints a and b. Before introducing p(°°), we perform one 
more transformation of the RH problem for S in the next subsection. 

4.5 Third transformation S ^-^ P 

The following transformation will modify the jumps in the vicinity of a and b: the jumps on Si 
and T,2 will be removed in Ua and Ub- As a drawback, a discontinuity will appear on dUa and 
dUb, but the jump matrices on these boundaries will be close to the identity matrix for large n. 
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Define 



Piz) := { 



Siz) 



for z G C\([/aU?7feUSs), 



1 1 



S(z)p(^)(z)-i^(n2/3/6(2))"3-3 ^ 1^ Ij e-(x-f)-3 for z G C/b \ Sg. 
Then P is constructed in such a way that it has jumps on a contour 



(4.40) 



Sp := (S5 \ (Ua U ?7fe)) U aC/„ U dUb U [a, (4.41) 

as shown in Figure m We define (n2/3/^(z))- I'^s ^nd {n"^/^ fa{z))-^''^ in such a way that 
they have branch cuts on [a,b] and they are positive on (6,00) and (—00, a) respectively. The 
jumps inside the disks on M \ [a,b] and the hps Si,S2 are equal to the identity matrix since 
S{z)P^^\z)-^ and 5(z)p('^)(z)-i are analytic there, but there is a jump on (a, b) due to the 
branch cuts of (?t-^^^/6(-z)) ^'^s g^j^j ^^2/3 j^^^-j-j-^o-a^ Also note that, unlike Y, T, S whose entries 
are all bounded in any bounded region of their domains, P{z) has inverse fourth root singularities 
at a and b. 

RH problem for P 

(a) P = {Pi,P2), where Pi is analytic in C \ Sp, and P2 is analytic in S'^ \ Sp, 

(b) we have 



where 



Jp{z) = < 



Js{z 
1 

V2 



P+{z) = P^{z)Jp{z), for zESp, (4.42) 

for ZGE5\(I4UI4), 



Le(x-f)-3 ( I ^1 ) (n2/3/^(^))i-3p(a)(^) for z G a;7a. 



^e(x-iV3 ( I (n2/3/,(z))l-3pW(^) fo, ^ e dUb, ^^'^^^ 



V2 



e-^ 
-e^ 



for z G (a, 6), 



(cl) as 2 ^ 00, Pi(z) =z'' + 0(2^=-!), 

(c2) as +00, P2 behaves as P2{z) = C'(e~('^+^)^), and as — )• 0, P2 behaves as P2{z) 

0{l). 

(c3) 



P(z) = {0{\z - a|-i/^), 0(|z - a\-^'^)) 
P{z) = mz-b\-'l'),0{\z-b\-'h) 



as z — )• a, 
as 2; — )• 6. 



(4.44) 
(4.45) 
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Figure 7: The contour Sp. On the boldface part of the contour, Jp = ^q"") and on the 
other parts, Jp ^ I uniformly. Note that Tip divides the complex plane into six regions: the 
two "edge regions" Ua and [7^, the two "bulk regions" in the upper and lower parts of the lens 
and not in Ua or Ub, and the two "outside regions". The dashed lines that belong to S5 but 
not to Ep, together with the interval {a,b), divide each edge region into four subregions, two 
inside the lens and two out of the lens. 

4.6 Construction of the outer parametrix 

For z G dUa U dUb, the definition of the local parametrices (|4.33p and ()4.38p together with the 
asymptotics (j4.30p for A imply that Jp{z) = I + 0{n~'^) as n — )• 00. For z G \ [a, h\ and not 
included in Ua or Ub, by the asymptotics of (t){z) given in ()4.23p and (|3.9p . we have that Jp{z) 
decays exponentially as n — )■ 00. Thus, in some sense, it is expected that 

P{z) (4.46) 

where ^'■'^•'(z) has the same analyticity, asymptotic, and periodicity properties, and has the 
jump condition 

)(x) = p1°°)(x)(_°^, 'q'^), forxG(a,6). (4.47) 
We would like to construct a solution to the following RH problem: 

RH problem for p(°°) 

(a) P('-) = {P['=°\p!f°^), where p[°°'^ is an analytic function in C \ [a, 6], and P^'^ is an 
analytic function in S"^ \ [a, 

(b) p(°°) satisfies the jump relation (j4.47p . 
(cl) as z ^ 00, P['^\z) = z^ + 0{z^-^), 

(c2) as — > +00, P2°^^ behaves as P2°°\z) = ©(e^^^^^^^), and as — > 0, P2°^^ behaves as 
Pt\z) = 0{l), 

(c3) 

P'^^\z) = {0{\z-a\~^/^),0{\z-a\-^/^)) as z ^ a, (4.48) 

pM(z) = (0(|^-6|-i/4),0(|z-6ri/^)) asz^b. (4.49) 

After the construction of p(°°), we will prove the convergence (|4.46p . 

We use the transformation Jci,co defined before in ()1.23p . where the parameters ci and cq 
depend on a and b, the endpoints of the support of the equilibrium measure. Recall Sa and Sb 
defined in ()1.24p and the relation p.4p between Sa,Sb and a,b. Below we write Jci,co &s J if 
there is no confusion. 
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By Proposition [21 J maps C \ D conformally to C \ [J(sa), J(sfe)], and maps D \ [— i, ^] 
conformally to S\ [3{sa),J{sb)], so that we can define the function F{s) on C\ (71 U72U [— ^, ^]) 

by 

)(oo)/ 



F{s) ' \ ^ (4.50) 

\P^)(J(.)) for.GZ)\[-i,i]. 

Since ^2°°'' is defined on §'^, that is, it satisfies a periodic boundary condition on S, we have 
that the definition of F{s) can be extended to (— ^, In this way the transformation from 
pioo) jr ig invertible: we can recover the outside parametrix by the formula 

Pi'^\z) = F{Ii{z)), for z eC\[a,b], (4.51) 

P^°°\z) = F{hiz)), for z G S\ [a,6], (4.52) 

where Ii and I2 are, as defined in ()1.26p and ()1.27p . the inverses of J mapping C \ [a, 6] to C\D 
and to D\ [— ^, ^] respectively. All information about the vector- valued function is now 

carried by the single complex-valued function F, which is discontinuous on 71U72 by definition. 

From condition |(c2)| of the RH problem for and the definition of J, it follows that F 

has a removable singularity at — ^, and that F has a zero of multiplicity A;-|-lat|iffc>0,a 
removable singularity if k = —1, and a pole of order —k — lifk<—l. The inverse fourth root 
singularities of at a, b are transformed into inverse square root singularities of F at Sa, Sb, 
because J'(s) has simple zeros at Sa and Sb- In order to compute the jump relation satisfied by 
F, note that 

s -+- i 

gJ(s) =gCis+co^2^ (4.53) 

and 

F+{s) = Pi^\j{s)), F^{s) = Pi^\j{s)), for . e 71, (4-54) 

F+(s) = pg)(J(s)), F^{s) = Pi'^\j{s)), forsG 72. (4.55) 

It is now straightforward to verify the following RH conditions for F. 



RH problem for F 

(a) F is analytic in C \ (71 U 72) if k > —1, and analytic in C \ (71 U 72 U {^}) if A; < —1, 

(b) for s G 7i U 72, we have 

s + i 

F+(s) =-e^i'^+^» ^F^{s), forsG7i, (4.56) 

F+{s) =e-'=i^-^« r^-(*)' for s G 72, (4.57) 



s+ 2 



(c) we have the asymptotic conditions 



F{s) = els" + 0{s''-^), ass ^00, (4.58) 

F{s) = 0{{s-^)'^'), ass^i (4.59) 

F{s) = 0{\s - Sa\-^) ass^Sa, (4.60) 

F{s) = 0{\s - Sb\~^) as s^Sb- (4.61) 
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One can explicitly construct a solution F to the above RH problem: 



F(s) = I ^ V 'b) ' (4_62) 



^(S-Sa)(s-S6) 



where y^(s — Sa)(s — Sfo) is taken to be continuous in C\7i and \/ {s — Sa){s — Sb) ~ s as s — )■ oo. 

Note that F{s) and the function Gk{z) defined in ()1.30p are related by (upon expressing Sa 
and Sb by 

G.(») = |?? (4.63) 
^ ^ \eJWF(s) if s G Z). ^ ^ 

4.7 The convergence of P — p(°°) 




Figure 8: The contour J ^(Sp). The boldface part consists of 71 and 72, the solid part is E' 
and the dashed part is S". 

We will now apply the same idea as in the construction of the outer parametrix to P, and 
want to transform the RH problem to the s-plane using the transformation z = J{s), in such 
a way that P = {Pi,P2) is transformed to a single complex- valued function J^. Therefore we 
define J" on C \ J^^(Sp) analogous to (|4.50p : 

/^i(J(^)) if5GC\5andJ(s)^Sp, 
\P2{J{s)) ifsGZ)\[-i,i] and J(s) ^Sp. 

The inverse of this transformation is given by 

Pi{z) = T{h{z)), forzGC\Sp, (4.65) 

P2{z) = HHz)), for zGS\ Sp. (4.66) 

The jump contour of J- will consist of the inverse image of Sp under J. We can decompose this 
jump contour J~^(Sp) into three different parts: 71 U 72, the part in D and the part in C \ D 
as follows, see Figure [H) 

J-i(Sp) = S'uS"u(7iU72), where S' = Ii(Sp \ [a, 6]), S" = l2(Sp \ [a, 6]). (4.67) 

Similar to F{s), the definition of T{s) can be extended to [— |, |) because of the periodicity of 
P2 and its behavior as ?R.z — )• —00. The RH problem for P, however, no longer transforms to a 
scalar RH problem for -P(s). For s G 71 U 72, we still have the scalar jump conditions 

T+{s) = - e-^(')-F_(s), for s G 71, (4.68) 

T+{s) = e-^^'^T-{s), for s G 72, (4.69) 
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but on the other parts of the jump contour, the jump conditions become non-local. Since 
J^±(s) = Pi^-|-(J(s)) for s £ T,' and J^±{s) = i-2,±(J('S)) for s G S", where the orientation for 
E' and S" is that inherited from the orientation on Sp through Ii and I2, the jump conditions 
()4.42p for P transform into 



-^4 
-^4 



(s) = Jp,n(J(s))^_(s) + Jp,2l(J(s))^_(l2(J(s))), 

is) = Jp,12(J(s))J-_((Il(J(s)))) + Jp,22(J(s))-^-(s) 



for s £ S', 
for s £ T," , 



(4.70) 
(4.71) 



where Jp is the 2x2 jump matrix defined in ()4.43p . In other words, the boundary value 
J^+(Ii(z)) depends not only on J^_(Ii(z)), but also on J^_(l2(z)), and vice versa for J^+(l2(z)). 
For this reason, we will call the jump relations ()4.70p - ()4.7ip "shifted" jump relations, and the 
RH problem for a shifted RH problem, following the terminology of [25j. The asymptotic 
conditions for are the same as the ones for F. By conditions (cl) - (c3) of the RH problem 
for P, we have analogous to ()4.58p -- ()4.6ip that 



k-l\ 



Hs) = 0{{s 

m = o{\s- 



1 



F{s) = 0{\s-s,\-'''), 



1/2 ^ 



as s 

as s 

as s 
as s 



00, 
1 

2' 

Sa ) 
Sb- 



(4.72) 

(4.73) 

(4.74) 
(4.75) 



Since F{s) 7^ for all s G C \ (71 U 72 U {|}), while at \ the order of the pole of T{s) is at most 
equal to that of F{s), we can define the analytic function 



F{s) 



for s G C\ J"^(Sp). 



(4.76) 



By (|4.68p . ()4.69p and ()4.56p . (|4.57p . it follows that R is analytic across (71 U72). Furthermore, 
the RH problem for F{s) and the shifted RH problem for -F(s) yield the following shifted RH 
conditions satisfied by R. 



Shifted RH problem for R 

(a) R is analytic in C \ (S' U S"), 

(b) R has the jump conditions 

R+{s) = Jp,ii(s)i?_(s) + jR,2i{s)R.ih{J{s))), for s G S', (4.77) 

R+{s) = Jp,i2(s)i?-(Ii(J(s))) + Jr,22{s)R-{s), for s G S", (4.78) 

where 

JrMs) = JpMHs)), JrMs) = JpMHs))^^^'^^^, (4.79) 

F{s) 

JrMs) = JpMHs))^^^^j^jy^, JrMs) = ^P,22(J(s)), (4.80) 

(c) R is bounded, and R{s) = 1 + 0{s^^) as s — )• 00. 
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Substituting the asymptotic properties of Jp stated in the beginning of Section 14.61 and the 
formula ()4.62p of F{s) into ()4.79p and ()4.80p . as n — t- oo, we have the uniform asymptotic 
estimates 

jR,ii{s) = l + 0{n-^), Jfl,2i(s) = O(n-i), forsGE', (4.81) 

jR,i2{s) = 0{n~^), Jij,22(s) = l + 0(n-i), forsGS". (4.82) 

Moreover, for s on the real parts of E' and S", J 11,21 vanishes identically: by ()4.43p and ()4.2ip . 
we have 

Jr,21 {s) = 0, for s G (S' U S") n M, (4.83) 

To obtain asymptotics for -R(s), we introduce an operator that acts on functions defined 
on T^R = S' U S". Let / be a complex-valued function defined on T,r. Then we define g = Ar/ 
by 

g{s) = [JR,n{s)-l]f{s) + JRMs)fiMHs))), forsGS', (4.84) 

5(s) = ^i?,i2(s)/(Ii(J(s))) + [J,j,22(s)-l]/(s), forsGS". (4.85) 

For bounded function /(s), g{s) is also bounded and decays rapidly as \s\ — >• 00. If we regard 
as a linear operator from L^(S/j) to itself, we will see that it is bounded and that its operator 
norm is 0{n~^) as n — )■ 00. For that purpose, note first that, by (|4.84p and (|4.85p . 

l|Aij/||L2(s^) < ||[J^,n - 1]/||l2(sO + II Jfl,2l/(l2(J))||L2(E') 

+ ll^i?,12/(Il(J))||L2(E") + Ilfe2 - 1]/||l2(s')- (4.86) 

Using the fact that Jr^u — 1 and Jij,22 — 1 are uniformly 0{n^^) on S' and S" as n — > 00, see 
(j4.8ip - (|4.82p . we obtain that there exists a constant c > such that 

W^RfhH^n) ^ ^WfhH^n) + \\JR,2lfiUm\LH^') + Pil,l2/(Il(J))||L2(E")- (4-87) 
For the second term on the right-hand side, we have 

||Jii,2i/(l2(J))|li2(2.) = [ U{h{J{s)))WRMs)fds 

= / |/(u)|Vi?,2i(Ii(J(u)))|2|(Ii(J))'(u)|dn (4.88) 

JS" 

< sup (|J^,2l(Il(J(n)))|2|(Ii(J))'(u)|)-||/||i.(^^). 

For u G Ti" bounded away from ±1/2, it is straightforward to verify by ()4.8ip and properties 
of the transformation J that | Jij^2i(Ii(J(^^)))Pl(Ii(J))'(^)| is 0(n~^) as n — )• 00, uniformly in 
u. For u G S" close to ±1/2, we observe by ()4.83p that Jij^2i(Ii(J('u))) = 0, which implies the 
existence of a constant c' such that 

\\jR,2if{hm\LH^n < ^II/IIl2(s,)- (4.89) 
Regarding the last term in (j4.87p . 

||JK,i2/(Ii(J))|li2(s.) = [ \f{h{3{s)))f\jRMs)fds 

JT." 

= I |/(n)|ViJ,i2(l2(J(n)))|2|(l2(J))'(n)|dn (4.90) 
Jt,' 

< sup (|J^,i2(l2(J(n)))n(l2(J))'(n)|) • ||/|li2(s,), 



28 



and it follows from (|4.82|) that 



||Jh,12/(Ii(J))IIl^(E") < -|I/IIl2(e^)- (4.91) 

From the above estimates, it follows that there exists a constant M > such that 

M M 
W^RfhHE,) < W^Rh^i^n) < -■ (4-92) 

Next, we define another bounded linear operator Ca^ from L^(S^) to itself, by 

Ca«(/):=C_(A^(/)), where g{s) = ^ lim [ /^d^, (4.93) 

and the limit s' — )• s_ is taken when approaching the contour from the minus side. The operator 
norm of Cajj is also uniformly 0{n~^) as n — >• oo since the Cauchy operator C_ is bounded. 
Thus (1 — Cajj) can be inverted by a Neumann series for n sufficiently large. We claim now 
that R satisfies the integral equation 

Ris) = l + C{AnR-)is), where Cg{s) = ^f j^d^ (4.94) 

To prove this claim, note that the solution to the RH problem for R is unique because it is 
equivalent to the uniquely solvable RH problem for Y. This means that it is sufficient to prove 
that the right-hand side of ()4.94p . which we will denote by R for simplicity, satisfies the RH 
conditions for R. Obviously R{z) is bounded and tends to 1 as z — )• oo, and it suffices to prove 
that the solution R satisfies the jump relations ()4.77p and ()4.78p . Using the Cauchy operator 
identity C+ — C_ = 1, it follows that 

R+-R- = {1 + C+{ArR^)) - (1 + C_(A^i?_)) = {C+ - C_)(ARi?_) = ArR^, (4.95) 

which implies indeed that R satisfies the jump relations (|4.77p and (|4.78p for R. Hence we 
conclude that R = R, and (|4.94p is proved. Since R satisfies (j4.94p . we have, taking the limit 
where s approaches the minus side of T^r, 

R--l = C-{ArR^) = Ca^R- - 1) + C_(A^(1)). (4.96) 

By the invertibility of (1 — Cajj), (|4.96p implies 

i?_ = l + (l-CA«)-iC_(A^,(l)). (4.97) 

This further implies that 

- MIlHSr) = Oin-'^), asn^oo. (4.98) 
Substituting (j4.97p into ()4.94p . we obtain an expression for R: 

R = l + C{ARil + (1 - Ca^)-1C_(A^,(1)))). (4.99) 
For s at a small distance 5 > away from the contour T,r, (|4.94p reads 

_^ir A^iR m) ,^ , 1 f Mm,,. (4,00) 
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The second term at the right-hand side of the above equation can be estimated by 0{6~^n~^), 
using the definition of the operator and asymptotic properties of Jr. Using in addition 
the Cauchy-Schwarz inequahty apphed on the first term on the right-hand side of the above 
equation, by ()4.98p we obtain 

\R{s) - 1| < ^\\^r{R- - 1)IIl2(s,) • II^IIl2(e«) + 0{6-'n-^) 

Although the estimate (|4.1U1|) does not work well for s in a 5- neighborhood of S/j, we note that 
for such s, given that 5 is small enough, the jump contour Tir can always be deformed in such 
a way that s lies at a distance 5 away from it. After this deformation, the above argument can 
be applied to obtain the uniform estimate 

R{s)-l = 0{n-^), as n cx), s G C\ Sh. (4.102) 

Through ^TjG^), (|i3T|) . and (fOSjl . the uniform estimate (|4.102!) yields 

Pi{z) = {l + 0{n-^))Pi°°\z), as oo, for z G C\Ep, (4.103) 
P2{z) = {l + 0{n-^))P^'^\z), as oo, for z e S'^VSp. (4.104) 

(n) 

The asymptotics for R can be used to obtain asymptotics for the polynomials pn by inverting 
the transformations 

Y^T^S^P^R. (4.105) 

We will do this in Section [6l 

5 Asymptotic analysis for the type I multiple orthogonal poly- 
nomials 

In a similar way as for the type II multiple orthogonal polynomials pj^\z), in this section we 

construct a RH problem for qj"\e^), and we analyze this RH problem asymptotically when 
j = n + k. Both the RH problem and the asymptotic analysis show many similarities with the 
ones for the type II polynomials, and once again the use of the transformation J will turn out 
to be crucial. 

In this section, we write qj^\x), the monic polynomials that define the multiple orthogonal 
polynomials of type I and satisfy the orthogonality relations ()1.12p . as qj{x) if there is no 
confusion. 

5.1 RH problem characterizing the polynomials 

Consider the Cauchy transform of qj(e^), 
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Due to the orthogonality (|1.12p . as z — )• oo, 



27riz J^l — s/z 



(5.2) 



= 0{z~^~^). 
For X G M, Cauchy's theorem imphes 

(5.3) 

Similar to (fi3]) . let 

X(z)=X(^'")(z):=(<7,(e^),Cg,(z)). (5.4) 
One verifies that X satisfies the following RH problem. 

RH problem for X 

(a) X = (Xi,X2), where X2 is an analytic function defined on C \ M and Xi is an analytic 
function on S'^, 

(b) X has continuous boundary values X± when approaching the real line from above and 
below, and we have 

X^{x) = X^{x) ^ j, forxGM, (5.5) 

(cl) as z — 00, X2 behaves as X2{z) = 0{z^^^^), 

(c2) as -)• 00 (i.e., Kz — )- +00), Xi behaves as Xi{z) = e^^ + 0{e^^~^^^); as — )- (i.e., 
5Rz — )• —00), Xi remains bounded. 

In an analogous way as for the RH problem for Y in Section HTTl it can be shown that X = X^-^'") 
given by ()5.4p is the unique solution to this RH problem. 

We will now perform an asymptotic analysis of the RH problem for X = as n — )• 00, 

with k a constant integer. This method will be to a large extent analogous to the nonlinear 
steepest descent method done in the previous section. Again we will construct a series of 
transformations of X and end up with a shifted small-norm RH problem. In order to emphasize 
the analogies with the previous section, we will use notations T, S,P, R, . . . for the counterparts 
of the functions T, S,P, R, . . . used before. 

5.2 First transformation X T 

Recall the functions g{z) and g{z) defined in ()1.33p . and define 

Tiz):=e-^Xiz)(^^ ^ eng(.)J e^^^- (5-6) 
Analogously to T in Section 14. 2[ T satisfies the RH problem 
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RH problem for T 

(a) T = (Ti, T2), where T2 is analytic on C \ M and Ti is defined and analytic in S'^ \ M, 

(b) T satisfies the jump relation 

f+{x) =f^{x)Jf{x), forxGM. (5.7) 

with 



Jf{x) 



gn{g_(x)-g+(x)) ^n{g-{x)+g+(x)-V{x)- 
{a;)-g- (a;)) 



(5i 



(cl) as z ^ 00, fa behaves as 12(2) = C'(z-(^+i)), 

(c2) as — > 00, fi behaves as fi{z) = e*^^ + ©(e'^'^"-'^)^), and as — )■ 0, fi behaves as 
fi{z) = 0{l). 

5.3 Second transformation T S 

By p.Sp . we have, like (|4.15p . the following factorization on [a,b]: 



Jf{x) 



1 



0\ / 1 



n4>-{x) ij \_l q) \^-n4>+(x) ^ 



1 



(5.9) 



where (p is defined in ()1.34p . Recall the lens defined in Section [4.31 and shown in Figure [5j 
Similarly as in ()4.19p for S, let us define S by 



Siz) := <^ 



Tiz) 
f{z) 



1 
1 0^ 



outside of the lens, 

in the lower part of the lens. 



(5.10) 



T(z) I _n,j,[z) ™ upper part of the lens, 



where 4){z) is defined in ()1.34p . Then similar to the RH conditions satisfied by S in Section [33 
we have the RH problem for S as follows. 

RH problem for 5 

(a) S = {Si, S2), where ^2 is analytic in C \ S5, and Si is analytic in S'^ \ S5, 

(b) we have 



where 



S^{z) = S-{z)J^{z), for z G S^, 
1 0\ 



(5.11) 



Jq{z) = < 



l\ 

-1 oj' 

.0 1 J 



for z G Si U S2, 
for z G (a, 6), 
for z e M \ [a,b], 



(5.12) 



(cl) as 2; — )• 00, §2 behaves as ^2(2) = ©(z-^'^+i)), 

(c2) as — )• 00, S*! behaves as S'i(z) = e'^^ + ©(e^*^"^)^), and as — )• 0, behaves as 
Si{z) = 0{l). 



32 



5.4 Construction of local parametrices near a and b 

In a similar way as for the construction of P^*^-* and P^^^ in Section [4.41 we can construct local 
parametrices -P^"^ and P^^^ in sufficiently small neighborhoods Ua and Ub of the endpoints a 
and b in such a way that they satisfy exactly the jump conditions 



Pi^\z) = Pl->{z)J^{z), 
Pf (.) = pf (z)Ja(z), 



p(a)/ 



z G S5 n Ua, 



(5.13) 
(5.14) 



Similar to the P^"-* (z) and P^^^ (z) defined in (j4.38p and (j4.33p respectively, the local parametrices 
P('^)(z) and P(^)(z) are expressed by 



P^-\z) := A(n2/3/a(^))e~^*^^)"^ 
pW(z) :=^(n2/3/,(z))e-t<^(^V3, 



(5.15) 
(5.16) 



where the functions fa and fb are as in (j4.32p and (|4.37p . A is as in (|4.26p . and the neighborhoods 
Ua and Ub as well as the contour can be taken the same as in ()4.38p and ()4.33p . We omit 
the details of the verification of ()5.13p and ()5.14p here, since almost identical arguments were 
used in Section 14. 4[ 



5.5 Third transformation S ^ P 

Define analogously to P{z) in (j4.40p . 

'S{z) 

S{z)P^-){zr^^^[n'l^fa{z))-\ 



P{z) = I 



for z G C \ {Ua U f/fe U S5) 



S(z)pW(z)-i^(n2/3/,(z))-i 



l]^' foi Z€Ua\^S, 

i'^^ ' ^ |j e^x'^a for z£Ub\ S5. 



Then like the RH conditions satisfied by P, P satisfies the following RH conditions. 



(5.17) 



RH problem for P 

(a) P = (Pi, P2), where P2 is analytic in C \ Sp, and Pi is analytic in S'^ \ Sp, 

(b) we have 

P+{z) = P^{z)Jp{z), forzeSp, 
where Sp is the same as in (j4.4ip . and 

'Jgiz) for ZGS5\(C^UI4), 



Jp{z) = < 



V2 



1 TTl „ 



1 

-1 0, 



1 1 

1 -1' 
1 1 



(n2/3/^(z))i'^3p('')(z) foTzedUb, 



for z G (a, 6), 



(5.18) 



(5.19) 



(cl) as z ^ 00, P2{z) = ©(z-^^'+i)). 
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(c2) as — )■ oo, Pi behaves as Pi{z) = e^^ + 0{e^^ ^)^), and as — )• 0, A behaves as 
Pi{z) = 0{l). 

(c3) 

P{z) = {0{\z-a\-^l^),0{\z-a\-^'^)) as z ^ a, (5.20) 

P{z) = {0{\z-b\-^/^),0{\z-b\-^/^)) as z^b. (5.21) 

5.6 Construction of the outer parametrix 

The RH problem for P has, as P, the property that its jump matrix tends to the identity matrix 
uniformly as n — )■ cx), except on [a,b]. We will first construct a solution to the following RH 
problem for p(°^) , which is the limiting RH problem (formally, ignoring small neighborhoods of 
a and b) for P as n — )• oo. 

RH problem for p(°°) 

(a) p(°°) = [P^°^\ P2°°^ ) , where P2°^^ is an analytic function in C \ [a, 6], and P^°°^ is an 
analytic function in S"^ \ [a, 5], 

(b) P(°°) satisfies the jump relation 

p|°^)(x) = p1°°)(x) (^_°^ J^, forxG(a,5), (5.22) 

(cl) as z ^ oo, Pi°^\z) = 0(z-('=+^)), 

(c2) as ^ +oo, P^°°^ behaves as p|°°^(z) = e^^^ + ©(e^'^"^)^), and as ^ 0, Pi"^^ behaves 
as Pi(~^(z) = 0{1). 

(c3) 

pM(z) = (0(|z-a|-i/^),0(|z-ar^/^)) as z ^ o, (5.23) 

pM(z) = (O(|z-6|-i/4)^0(|z-6|-^/^)) asz^6. (5.24) 

Inspired by the construction of p(°°) in Section IT6l we search for p(°°) in the form p(°°) (z) = 
{F{l2{z)), F{Ii{z))), where Ii and I2 are, as before, the two inverses of the map J defined in 
(frkjl and ^rm . Hence 

fpi°°^(J(s)) for s G C\P», , , 

[P| '{J{s)) for s G P»\ 

and like P(s) in Section WM F can be analytically continued to [— ^, ^). At ^, F has a pole of 
order k ii k > 0, a removable singularity if A: = and a zero of multiplicity —k if A; < 0. From 
the RH conditions for p(°°), we deduce the following RH problem for F. 

RH problem for F 

(a) F is analytic in C \ (71 U 72) if /c < 0, and analytic in C \ (71 U 72 U {^}) if A; > 0, 

(b) for s G 7i U 72, we have 

P+(s) = P_(s), forsG7i, (5.26) 

F+is)= -F_(s), forsG 72, (5.27) 
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(c) we have the asymptotic conditions 

= 0(s-('=+^)), ass^oo, (5.28) 

F{s) = e^(i'=i+'=o)(. - 1)-'= + 0{{s - as s ^ ^, (5.29) 

F{s) = 0{\s - Sa\-^), as S^Sa, (5.30) 

F{s) = 0{\s - Sb\~^), as s^Sb. (5.31) 

It is verified directly that 



F{s) = ^ i ^/ ^ e^(2^^+^o)(g - )-\ forsGC\72, (5.32) 



solves the above RH problem. Here we choose the branch of the square root y^(z — Sa){z — Sb) 
that is analytic except on 72 and close to 2; as 2; — )■ 00, 

5.7 The convergence of P — t- 

Define analogous to J^{s) in (|4.64p 

^(.) := if..C\I.andJ(.)^Sp, 

\Pi(J(s)) ifseP\[-i,i] andJ(s)^Sp. 

We have the scalar jump conditions 

J"+(s) = J"_(s), forsG7i, (5.34) 

T+{s) = - /_(s), for s G 72, (5.35) 

and the shifted jump conditions 

f+{s) = Jp ,i(J(s))A(s) + Jp 2i(J(s))A(Ii(J(s))), for s G S", (5.36) 

T+{s) = Jp^^^{3{s))f^{{h{J{s)))) + Jp^^^{3{s))T.{s), for s G S'. (5.37) 

The asymptotic conditions are the same as those for F{s) 

f{s) = ©(s-^'^+i)), as s ^ 00, (5.38) 

^(s) = e^(^^i+^»)(s-i)-'= + 0((s-^)-^+^), ass^^, (5.39) 

:F{s) = 0{\s-Sa\-^^^), ass^Sa, (5.40) 

J"(s) = Ods - s^r^/^), as s^Sfe. (5.41) 

Next we define, analogous to R{s) in (j4.76p . 



R{s) := j^, for s G C \ J"^(Sp). (5.42) 



Then like R, R is analytic at ^ and across (71 U 72), and satisfies the following shifted RH 
problem. 
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Shifted RH problem for R 

(a) R is analytic in C \ (S' U S"), where S' and S" are defined in (j4.67p . 

(b) R has the jump conditions 

R+{s) = J^^^^{s)R^{s) + J^_2i(*)^-(Ii(J(^)))' foi- ^ e S"' (5-43) 

= J^^,2(«)^-((l2(J(s)))) + Jji^22i^)RAs), for s E S'. (5.44) 

where 

= ^Pii(J(^)), JR2ii') = Jp2ii^i') f^^f['^^^ ^ (5-45) 

• ' ' ' F{s) 

jR,ui^) = Jp,l2i^(^) f^^f['^^\ Jr,22{^) = Jp^^i^))- (5-46) 

F{s) 

(c) R is bounded, and R{s) = 1 + ©(s"^) as s — > oo. 

As n — 7- c«, we have the uniform asymptotic estimates analogous to ()4.8ip and ()4.82p 

J^ ii(s) = l + J^2^(s) = O(n-i), for s G S", (5.47) 

Jr,12^^) = JR.22i^) = 1 + C'(ri"^), for s G S', (5.48) 

These estimates imply, in a similar way as ()4.8ip and (j4.82p do in Section 14.71 the uniform 
convergence of i? to 1: 

R{s) = l + 0{n-^), as oo, for s G C\ (S'US"). (5.49) 

Hence, by ([02]) . ([OH]) , and (lOHp . we have, like (|TO3]l and (ICTHp . 

Pi(z) = (1 + C'(n-i))Pi^°°^(z), asn^oo, for zGS'\Sp, (5.50) 
P2(z) = (1 + C'(n-i))P2^°^^(z), asn^oo, for zGC\Sp. (5.51) 

6 Proof of main results 

In this section we collect the asymptotics of p^^l^j^{z) and q~^^f^{e^), from the analysis in Sections 
[Hand El The goal is to prove Theorem [21 

6.1 The asymptotics of p'^]^f.{z) 

(n) 

The main task in the computation of the asymptotics for p^^^f, consists of the inversion of the 
transformations F H> T H> 5 H> P. By ([i3]l . ([4.13p . (|4.19p . ([4.40p and the asymptotics obtained 
in Section [4.61 we will find the asymptotics of p^"^^. In Figure [7] it is shown that the complex 
plane is divided into the outside region, upper and lower bulk regions and two edge regions by 
Sp. We restrict ourselves to the upper half plane because of symmetry, and do the computation 
in each of the four regions. 
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Outside region For z in the outside region, we have 

Vn\{z) = y/^+'^'^^lz) = ri(z)e"g(^) = Si(z)e"g(^) = Pi(z)e"g(^). (6.1) 
By (jTO3]l and (li3TD . 

pSfc(^) = (1 + 0(n-i))F(Ii(z))e"g(^), as n ^ oo, (6.2) 
where F is defined in (|4.50p . Substituting the identity (j4.63p for F into (|6.2p . we have 

pSfe(^) = (l + 0(n-i))Gfc(Ii(z))e"sW, as n ^ oo. (6.3) 
This proves (|1.37p for z in the outside region. 

Bulk region Similar to (|6.ip ~ (|6.3p . for z in the upper part of the lens but not in Ua and ?7ft, 
we obtain 



"g(z) 



= {Si{z) + 52(^)e~"'^(")+^)e"sW 

= Pi(z)e"s(^) + P2(2)e^e"(^(^)-s(^)+^) 

= (1 + 0(n-i))F(Ii(z))e"g(^) (6.4) 
+ (1 + 0(n-i))F(l2(z))e^e"(^(^)-s(^)+'^) 

= (l + 0(n-i))Gfc(Ii(z))e"g(^) 

+ (1 + 0(n-i))Gfc(l2(2))e"(^(^)-^(^)+^), 

as n — ;> oo. In the last identity of ()6.4p we use (|4.63p and the identity = e'^i^2(2)+co Mfl±|.^ 

which is a consequence of ()1.23p and the definition of I2 as an inverse of J. We obtain the 
formula (jl.39p for z in the upper bulk region. 

In particular, if x G (a, b) and z — >• x from above, we have by ()3.8p that V{x) — g+(x) + £ = 
g_(x), and further from the definition (jl.33p of ^{z), we have g±(x) = Jlog|x — y\d^v{y) i 
TTz J^dfiv On the other hand, as z — )• x from above, by ()1.28p and ()1.29p . 11(2;) and 12(2) 
converge to I+(x) and I„(x) respectively. Noting that I_(x) = I+(x), we have from (j6.4p and 



p2fc(x) = rfc(x)e"/i°sl^-^l'^^^(f) cos (^n^ j\i^v{t) + 0fc(x)) + 0{n-^] 



(6.5) 



where rfc(x) and 9k{x), as defined in ()1.3ip . are the modulus and argument of 2Ga,.(I+(x)) 

^^fc (I+(x)+i)(I+(x)-i)^ 



'V(I+W-^a)(I+(x')-S6)' 

Edge region For brevity we only consider the case z E Ub, the case z € Ua can be treated 
similarly. As shown in Figure [TJ the part of Uh in the upper half plane is divided by the lens 
S5 into two parts, one in the lens and one out of the lens. If z G UbCi C"*" is outside the lens, 
we obtain 

p(j^l^{z) = y/"+'=''^)(z) = Ti(z)e"s(^) = 5i(z)e-^(^), (6.6) 
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and by (Onj) . 



1 1 
-1 1 



6 4 2n6/b(z)4 -6 4 2n 6/f,(2:) 

e 4, ^ 2 jh{z) i e 4^2n e fi,{z) 4 



11 ^(^^(^). 



(6.7) 



By (I03]l . and (lilIl3|) - (liTi1i]l . we further obtain 



(Pi(z) - ^P2(^)e") n6/,(^)4 Ai(n3/,(z)) 
(Pi(z) + iPaWe^) n-hbiz)--^ Ai'(ni/b(z)) 
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(g{2)-g(2)+y(^)+^) 



^(Ii(z)) - i^(l2(z))e-^(^^(^))) A(z)3 Ai(nf 



n-6/,(z)-4 Ai'(n3/,(z)) 



-F(Ii(z))+i^(l2(z))6-^(^2(")) 
(gfc(Ii(z)) - igfc(l2(^))) n^/,(z)t Ai(nt 
(gfc(Ii(z)) + igfc(l2(^))) n-^/b(^)-3 Ai'(ni/,(z)) 



gf(g(2)-g(^)+y(^)+^) 
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(g(2)-g(2)+y(^)+£) 



where Gk is defined, analogous to the formula ()4.63p for Gk, as 

C / ^ ._ /•^(*) if s € CVL) and J(s) ^ Sp, 

- |eJW^(5) if s G D \ [-1, 1] and 3{s) ^ Sp. 

From (|I775|) and (|4.1U2p . we have that 



as n — )• oo. 



(6.8) 



(6.9) 



(6.10) 



Hence we obtain ()1.47p for z in the edge region Ub, upper half plane, and outside of the lens. 

Let us now focus on the asymptotics of p^j^^i^{z) for z = 6 + f'^{b)~^n~'^/^t which is in the 
upper half plane and outside of the lens, where t is bounded. Then 

Ai(n2/3/6(z)) = Ai(t) + 0(n-2/3)^ 

as n — 7- oo, 

Direct computation yields, as n ^ oo, by (fOHD - ffTMl) and (fOejl - if]^ . 



lAz) 



Sb + 



Sb 



hiz) = Sb 



(gb + ^)(g&-^) ^;^^^_l/2 



n 



Sb 



and that as s — )• Sb, by ()6.10p and (|1.30p . 
g,{s) = (1 + O(n-i)) (2-H] + -r'^ct" 



1 



1, 



4 ^ ci 2 



t(l + 0(n-2/3)f), 
t(l + 0(n-2/3)t), 



+ 0{s-sb) 



k-l 



'4 61 

where all square roots take the principal value. Hence when t is bounded 

'1 (\ 



1 



\/s- Sb 



(6.11) 

(6.12) 
(6.13) 

(6.14) 



Qk{\^{z))^iGk{\2{.z)) = 0{n-\)t\. 



1 1 1 

4 ^ ^ ~ 2 



fc-i 



c 



+ 0(n- 



1 _i 
net 4 . 



(6.15) 
(6.16) 
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Substituting (j6.15p and (|6.16|) into (j6.8|) and noting that A (6) = and /^(&) > 0, we obtain 
()1.49p for z outside of the lens. 

If z G f/fc n C"*" and inside the lens, then like 



Vnlk^z) = {Si{z) + S2(z)e-"<^(^)+^)e"g(^), (6.17) 

and like (fSTTll . 

(5i + S2e-"<^(^)+^ ^2) = 

^/2(A,P2)e(x-f)-(_l^ |)''(n^/V.(.))^-P(^)(^)(^_„i)+. °) • (6.18) 

Hence by ([051) - (|I3T|) and (|4.103p ~ (|4.1U4|) . and using the identity A\{x)+uj Ai{ujx)+uj'^ Ai{ij?x) = 
0, we find that the result in (|6.8p still holds, and so do the subsequent asymptotic formulas 
I^M)-<^JB- Thus we can still prove (fTITll and (fOUD . 

6.2 The asymptotics of q:^l^{e^) 

The derivation of the asymptotics for ^Z^+fc(e^) is similar, and we need to invert the transfor- 
mations X I— r I— 5 I— )■ P using ()5.4p , (|5.6p , ()5.10p , and ()5.17p . For brevity, we only consider 
the outside region and the bulk region. 

Outside region If z is in the upper half plane and not in the lens or Ua, Uh, we have 

g(^^^(e^) = +^'")(z) = fi(z)e"s(^) = Si(z)e"^(") = Piiz)e''^^'\ (6.19) 
By (f03]) and ([532]), we find similar to (^J^ 

q^:Uen = ^(l2(z))e"^(^) = R{h{z))F(h{z))e-^(^\ (6.20) 
By the formula (|5.32p for F and the asymptotic formula (|5.49p for R, this yields 



Ckiel = (1 + 0(n-^))^=^^^^=e'=(l-+-)(l2(.) - V^^e"^(^), as n ^ oo. 

(6.21) 

In ()6.2ip and later in ()6.23p . \J {z — Sa){z — Sb) is chosen to be close to z as z — oo and has 
branch cut along 72. Substituting Sa and by (jl.24p . we prove (jl.38p for z in the outside 
region. 



Bulk region Similar to ()6.4p . 

gj^^^(e^) = fi(z)e"^(^) = 5i(z)e"sW + 52(z)e"(^(^)-s(^)+^) = A(2)e"^(^) +P2(2)e"(^(')"S(")+^). 

(6.22) 
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By ([533]), (1^^ and (fOUD . we find that as n oo, 

= i?(l2(z))F(l2(^))e"^(-) + i?(Ii(z))F(Ii(z))e"(^(^)-s(^)+^) 

= (1 + 0(n-i))^^^^^^^=e'=(5-i+^«)(l2(z) - i)-'=e"g(^) 



+ (l + 0(n-i)) 



a/^^ "^^^ .fc(|^i+^o)fTir.) - l)-fcen(y(.)-g(.)+£) 

V(Ii(^)-Sa)(IiW-56) ^ ^ ^ 2^ 

(6.23) 



Substituting Sa and Sb by ()1.24p . we prove p.40p for z in the upper bulk region. 

As z — 7- X G M from above, noting that V{x) — g-(x) + i = g+{x) by (|3.8|) . 12(2) I- (2^)) 
and li{z) — )• I+(x), and using the identities = I+(x) and g±(x) = / log|e^ — e^|d;Uy(y) it 



d/xy, we have like (j6.5|) . 



cosfnvr / dfj,v{t) + 0k{x)] + 0{n ^) 



3.24) 



where rfc(x) and 9k{x), as defined in ()1.32p . are the modulus and argument of 2Gfc(I_(x)) = 
6.3 Proof of Theorem [1 

The asymptotic results obtained in the last two subsections nearly prove items |(a)[ |(b)| and 
part of (c) and 1(d) I of Theorem [2j However, in the statement of the theorem, the regions where 
asymptotic formulas are given, are Ag, Bs, Cg, and Ds, which are similar but not exactly equal 
to the outside, upper bulk, left edge and right edge regions that depend on Sp. We observe 
that if is a fixed small enough number, we can take the radius of Ua and Ub large enough so 
that they cover Cs and Dg. On the other hand, we can also take the radius of Ua and Ub small 
enough and the shape of the lens thick enough to let the upper bulk region cover Bs, and we 
can take the radius of Ua and Ub small enough and the shape of the lens thin enough to let the 
outside region cover As. In this way, by using different contours Sp, the asymptotic results in 
the outside, upper bulk, left edge and right edge regions are translated into results in regions 
As, Bs, Cs, and Ds respectively. 

Although we have not proved all the asymptotic formulas in items (c) and |(d)| of Theorem 
[21 the remainders can be proved using the method presented in the previous two subsections, 
and we omit the details. 

To compute h^^^j^ and prove item (e) of Theorem [21 we note that it appears in the leading 



coefficient of Cpn+k{z), see ((331) • Using ([32]), KWf . Kl9\\ and ([330]), we have, for z in S, 
outside of the lens and away from a and b, 



Cpn+k{z) = 



(n+fc,n) 



(^) 



-ng(-)r2(^ 



e"^e-"sWp2(z) 



Since ^{z) = z + 0{e~^) as z ^ +00 in S, (f05]) yields 

= _j^^-ni^-ik+i)z ^ o{e~^''+^>), as z ^ +00. 
2m 



(6.25) 



(6.26) 
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By dSSU, (HISS]), and ([322]), we have as n oo, 

lim P2(^)e(^'+^)^ = lim P^°^\z)e'^''+'^> {1 + 0{n-'^)) 
= lim F(l2(z))e('^+i)^'(l + 0(n-i)) 

/ 1 \ '^^^ -Cll2(2)-C0 

= Hm cMl2(z)-- -==_===e('=+^)^(l + 0(n-i)). 

^^+°° V 27 ^(l2(z) - Sa){h{z) - Sb) 

(6.27) 

From the formula ()1.23p of J(s) which is the inverse function of l2(-2), we have 

l2(z) = i + e^+'=Oe-^ + 0(e-2^), as z ^ +oo, (6.28) 

and we obtain that 

^hrn^c^ (l2(z) - i j ^ .(^+^)- = c^e'=(^+'=°) 



^{I2iz)-Sa){l2iz)-Sb) ^(i -Sa)(sb- i) 



(6.29) 



where Sq and are expressed in ci by (|1.24|) . Formulas (|6.29|) , (j6.27|) and (|6.26|) yield Theorem 

El 



A Proofs of several technical results 
A.l Proof of Proposition [T] 

Our proof is similar to that of [151 Proposition 2.1]. By the formula of the probability density 
function ()1.9p . the average of nj'=i(^^ ~ ^^^) *^an be expressed as 

n ^ „ n n 

,=1 ^nJR-j^^ ^.^1 

From this formula, it is clear that IE^(nj=i(^^ ~ ^^^)) a linear combination of e^^ with 
k = 0, 1, . . . ,n, and that the coefficient of e"^ is equal to 1 since the probability measure is 
normalized. To show that it is equal to gi"^(e^), we only need to verify that it satisfies the 
orthogonality conditions ()1.12p . which characterize q^\e^) uniquely. 

Expanding the Vandermonde determinant over the symmetric group Sfi gives 

n 

n(A, - A,) = det(Ar').,,=i,...,n = E ^-^y n 4'^''- (^.2) 
Substituting ([Q) into (EOl) . we obtain 

n ^ „ n n 

j=i ^"-^^".es,. i=i i<j 3=1 

, r. n n 



(A.3) 
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Substituting the identity 



Y[{e' - JJ(e^^ - e^') = det 

j=l i<j 



Ai 

1 



/I e 



1 e^" 



(A.4) 



into ()A.3p . we obtain after integrating with respect to Aj that 



/ "^00 



3 = 1 



det 



moi 
mil 



min 



\ 1 



_n2 



where 



(A.5) 



Then it is straightforward to verify that for k = 0,1,... ,n — 1, 



/ E;(n(e^-e'0)^V"^^^)dz = ^det 



( moo 

mio 



moi 
mil 



mn-l,0 "T-n-1,1 



nion \ 
min 

rrikn I 



0. (A.6) 



Thus we prove that IE^(nj=i(^^ ~ ^'^■' )) satisfies the orthogonahty condition ()1.12p that deter- 
mines ctn^(e^), and then it follows that IE'„(nj=i(6^ ~ ^^^)) = 9n"^(e^). 

A. 2 Proof of Proposition [2] 

In this proof, we fix ci G and cq G K, and J stands for Jci,co such that J(s) = cis + cq — 



log — f • Recall that s„ 
(ET 



""^ + ^ as in ()1.24p . and a = 3{sa), b = 3{sb) as in 



Cl 



To prove part (a) , we show that the equation J(s) = x: 

(1) has a unique solution s in the upper half plane C"*" = {s = u + iv with v > 0} if x € (a, 6), 

(2) has no solution in if x G R \ (a, 6). 

Moreover, as x runs from a to b, the solutions s = s{x) form an arc in from Sa to Sb- Then 
this arc is the desired 71 in Proposition [21 and the complex conjugate of 71 is the arc 72. 
For s = u + iv with > 0, J(,s) G M if and only if the identity for its imaginary part 



civ — arccot = 

V 



(A.7) 



is satisfied, where the range of arccot is (0,7r). It is a direct consequence of (|A.7p that v < ^. 
Under the condition < v < ^, (|A.7p is equivalent to 



u 



1 . N 2 

-T + V cot(cif ) — V 



(A. 
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By direct calculation we find that the right-hand side of (jA.Sp is a decreasing function in v for 
< f < ^. Moreover, as v ^ 0, it tends to | + ^, and as v ^ it tends to — oo. 

Thus for J(s) to be real where s = u + iv with v > 0, u has to be in (s^, Sfe), and for any u 
in this interval there is a unique v > to make (jA.Sp hold. The locus of all such s = u + iv is 
an arc in connecting Sa and s;,. As a consequence of (|A.8p . v increases as u runs from Sa to 
0, and then decreases as u runs from to Sb- At any s in this arc, 

and it follows that J is a homeomorphism from this arc to the interval [a, b], which proves part 



(a) of Proposition [2j 

Next we prove part |(b)[ It is easy to check that J maps the ray (s^, cxd) to (6, oo) and the ray 
(— oo, Sa) to (—CO, a) homeomorphically. Then it suffices to show that J is a univalent map from 
C+ \ D onto C"*", and the univalent property of J on C~ \ D follows by complex conjugation. 
To this end, we use the following elementary lemma: 

Lemma 3 (Exercise 10 in Section 14.5 of [l9]). Suppose that G and Q are simply connected 
Jordan regions and f is a continuous function on the closure of G such that f is analytic on 
G and /(G) C Jl. // / maps dG homeomorphically onto d^l, then f is univalent on G and 

f{G) = n. 

But this lemma is not directly applicable, since both \ D and are unbounded. Let 
di^) ■— be the conformal map from the unit disk to the upper half plane, we find that 
g^^ o J o ^ is a map from the simply connected region g~^{C~^ \ D) into the unit disk, and the 
map is homeomorphic on the boundary. A direct application of Lemma [3] shows that g^^ o3 o g 
is univalent in (^~^(C^ \ D) and onto the unit disk, hence J is univalent in \ D and onto the 
upper half plane, and part |(b)| is proved. 



To prove part (c) we find by direct calculation that J maps homeomorphically 

(1) the interval (sa, —\) to the ray (— oo,a), 

(2) the interval (|,Sf,) to the ray (6, oo), 

(3) the upper side of the interval {—\, \) to the horizontal line M — vri, and 

(4) the lower side of the interval (— ^, \) to the horizontal line M + vri. 

Then it suffices to show that J maps D n onto S n univalently. We use Lemma [3] again. 
Similar to the conformal map g^ we use the conformal map /i(s) := log |^ to transform the 
unbounded region SnC^ into a bounded one. We omit the details since the arguments are very 
similar to those in the proof of part |(b)[ 

A. 3 Proof of Lemma [2] 

First, we show that for any xi G ^ ()3.3p has a unique solution as an equation in xq. Note 
that 

^ ''V (J.,.o(^)) -^ = hh" i^^uM) ^ = ^9 / V" (J.,.o(.)) ^ 



2 •'T 2 •'Ti ^2 

- / ^ (j^i,^o(s(^)))9 — de, 

(A.IO) 
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where we parametrize s E 71 by its argument 6 that runs from to vr. This parametrization 
is weU defined since as s moves along 71, its imaginary part increases as its real part increases 
from Sa to 0, and then decreases as its real part continues to increase from to Sh, as shown in 
the proof of Proposition [2j 

Below we show that the right-hand side of (|A.10p is bounded below by a positive constant 
for all xq G M. Since ^"(Jxi,2:o('S(^))) is bounded below by a positive number by the strong 
convexity of V, we need only to prove for all 9 G (0, vr), 3'log(s(6') — ^) = aig{s{9) — ^) is 
an increasing function. We show the increasing for 6 € (f jTI") and 9 G (0, ^) separately. For 
geometric reasons, when 9 G (f ,vr), arg(s(^) — ^) is increasing with 9 since both ^s{9) < and 
^s{9) > are decreasing. For 9 G (0, we use the identity 

91og(s(0)-i) = 9(^xis(^) + :Eo + log(s(^) + i)) -9J,„,Js(^)). (A.ll) 

Here QJxi,xo{s{9)), by the construction of 71, vanishes, Qs{9) increases as 9 runs from to ^ 
and for geometric reasons 9 log(s + ^) also increases as 9 runs from to ^. Thus we have that 
for 9 G (0, f ), '=slog{s{9) - |) = xi3s(6') + 91og(s(6l) + i) is increasing. 

Now we have that as a function in xq, ^ !f^^' {'^xi,xo{s)) / {s — ^)ds is a bijection from M 
to M, since its derivative is bounded below by a positive constant. Hence by continuity, there 
must be a unique xq to make this function equal to 1. Given xi G M"*", we denote the unique xq 
that solves ()3.3p by co{xi). Similarly we can show that co(xi) is a continuous function in xi. 

Although we do not have a simple formula for cq{xi), we show below that 

- — : ^ V' {J x^^co{xi){s))ds < x^^, for xi sufficiently small, (A.12) 

^ j) V'{Jx^^co{xi){s))ds > for xi sufficiently large. (A.13) 

Hence, by continuity, it follows that there exists ci G that, together with cq = co(ci), solves 

As xi — 7- 0^, from (|A.8p . it follows that the shape of 7 is close to the circle with radius x-^ 
and center 0. Hence if we parametrize s G 7 as before by its argument 9, we have for 9 G [0, 27r), 

s{9) = e^'x^'^' + 0(1), hm = i, lim V {3x,,xo{s{9))) = V'{xo). (A.14) 

xi~^o+ s{9) — ^ 2:1^0+ 

By the dominated convergence theorem, we have 

r2TT 



lim — j> V'{Jx,,x,{s)) r = V- / y'iJx.,xo{s{9)))—^ds 



1 

2Tri 



2n 



(A.15) 



V'{xo)ids = V'{xq). 



We find lim^^_^o+ ^0(2:1) = xq, where xq is the unique value such that V'{xo) = 1. From the 
results obtained above, we have that 

77- / V'{Jx,,co{xi){s))ds = o{x^^) as xi 0+, (A.16) 
Ztti 

— 1/2 

since the shape of contour 7 approaches to the circle with radius x-^ , and the integrand tends 
uniformly to V'{xq) = 1. 
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On the other hand, for large values of xi, we use the expression 
7^ lv'{3^,,^,{s))ds = — 9 / V'{i^,^,,{s))ds = — r V'{3^,,^Mu)))^v'{u)du, (A.17) 

where s S 71 is expressed as a function in its real part u = "Sis, and v{u) > is defined by the 
condition that s{u) = u + iv{u) G 71, and Sq, Sh are the two endpoints of 71, as denoted in the 
beginning of Appendix lA.21 with the parameters ci,co substituted by xi,xo. Let us decompose 
the integral at the right of (|A.17p as Ii + 12 + /s, where 

h = — I ' V'{3^,,,„{s{^)))'^v'{u)du, (A.18) 

1 

h = — V'{J,,,,,isiu)))'^v'iu)du, (A.19) 
l3 = — r V'{3^,,,Mu)))'^v'{u)du. (A.20) 



From (jA.Sp . it is not difficult to find that as xi 



00, 



vi-\) = vi\) = l^i' + o{x^'). (A.21) 

We know that V' is an increasing function in u and that v{u) is an even function. From 
Appendix IA.2I we have that v{u) is increasing for u G (sa,0) and decreasing for u G (0,Sb). 
Hence the integral I2 is positive. Using the monotonicity of V' and integration by parts for Ii 
and I3, we similarly obtain 

h+h> ^V'i3,,,,Mh))ivih - vi^b)) - -V'i3.,,.M-h))ivi-h - v{sa)) 

^ ^ ^ 2 2 ^^_22) 



where in the last line we used the identities v{sa) = v{sh) = and v{—^) = v{^). Hence (|A.17p 
and the estimates of I2 and /i + I3 above imply that 

^.jv\3^,,^,{s))ds > i (v'iJ,,,M\)) - V'{J,„,m\)?) (A.23) 

As xi oc, 

1 X\ 1 X\ 

Jx^,xoisi--)) = xo- Y + oixi), Jxuxo{s{-)) = X0 + Y +o{xi), (A.24) 

where the two o(xi) terms are independent to xq. By (|A.24|) and the assumption V"{x) > c > 
for all X, we have that if xi is large enough, then uniformly for all xq G M 

V'iJx^xoisil))) - V^'(J.„x„(s(-^)) > CXI. (A.25) 
Substituting (Ia:25]) and into (1X23]) . we have that as xi — t- 00 and xq = co(xi), 

^ j( F'(J,,,,„(,,)(s))ds » xr^ (A.26) 

We note that f^V' {J^^^co(xi)is))ds is continuous in xi, since f^V {Jx^^xo{s))ds is con- 
tinuous in xi,xo and co(x) is continuous. Then we find that the estimates ()A.16p and (]A.26p 
imply that there is a pair (ci,co = co(ci)) such that both (j3.3p and (|3.2p are satisfied. 
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B Explicit construction of the equilibrium measure for quadratic 
and quart ic V 

In this appendix we use the method developed in Section[3]to find the endpoints of the support of 
the equiUbrium measure exphcitly for quadratic and quartic external fields V. In the quadratic 

2 

case, we consider a monomial external field V{x) = but the same method can be applied 
to all quadratic V. We also construct the density function of the equilibrium measure. In the 
quartic case, we confine our attention to V such that V{x) — | is an even function. Under this 
condition the equilibrium measure is symmetric around the origin. In contrast to the quadratic 
V that is automatically convex, we also consider quartic V that is one-cut but not convex. 

2 

External field V{x) = ^ 

In this case, V'{x) = j, and a simple calculation of residue yields 



(f) V' { ClS + Co - log 



s + 



ds = 

t' 2m 



V [cis + Co - log 



t 2t 



Thus by Lemma [21 we have 

t 

Ci =t, Co = -. 

The support of the equilibrium measure, as expressed by ()3.4p . is 

[a,b] = [Jf i(Sa), Jj i(Sa)] 



(B.l) 
(B.2) 



-{t - Vt^ + 4t) - log , -(t + Vt^ + 4t) - log 



B.3) 



In particular, for t = 1, we have 

. ^, r-\/5 + l , 3 + ^/5^/5 + 1 ^ 3-^/5' 
[a, b] = ^—^ log ^— , ^— - log ^— 

To find the equilibrium density, we have as a particular case of p.2ip that 



(B.4) 



M{s) 



^logy, forsGC\A 



for s £ D. 



Then by (|3.23|) . after a straightforward calculation, we obtain the following expression 

vr 



(B.5) 



(B.6) 



where 1+ is as before the boundary value of the inverse of J = i which parametrizes the 
curve 7i. 



46 



External field V{x) = x^/i + ux'^/2 + x/2 

In this case, V'{x) = x^ + ux + ^, and the calculation of residues yields 

^ / y ( CIS + Co - log ) ds = ^ + 3ci + 3cl + u, (B.7) 



1 cis + co-log^ ^3 3 3 3,2 ^ 3 1 

-ds = ^ + i-^ + -)c? + (-^ + 6co + -)ci + eg + nco + 



27ri/^ s-i 8 ■ ^ 4 ' 2''" ' ' 2 ' "'^ ' 2''^ ' ' ' 2' 

(B.8) 

As a consequence of the relation V{x) = V{—x) + x, the equilibrium measure //y is symmetric 
around the origin. Indeed, changing variables s i— )• — s and t i— )• — t in the energy functional 
(jl.l7p . it is straightforward to verify that lyifJ-v) = Ivii^'v)-, where Jly is defined by the fact 
that = ^y{—A) for any Borel set A. From the uniqueness of the equilibrium measure, it 

follows that /iy = /iy. In particular this implies that the support of the equilibrium measure is 
of the form [—6, b]. By (|1.21|) . we have cq = 0. Substituting this and ()B.7p into (|3.2|) . we obtain 
the equation 

cf + 12c? + 4uci -4 = 0. (B.9) 

Remark 8. Although the equilibrium measure, which is the limiting mean eigenvalue distribution 
of the random matrix ensemble as n — )• oo, is symmetric around the origin, this is not true for 
the finite n joint probability distribution of eigenvalues ()1.6p . The latter would only be invariant 
under the change of variables Aj — t- — Aj if the term x/2 in V{x) were replaced by — ^)x. 

For any value of u, the equation (|B.9p has a unique positive solution by Descartes' rule of 
signs. We have an explicit formula for ci G M"*" in u by the formula for the roots of a cubic 
equation, but we will not write down the long formula. Together with cq = 0, ci > gives us 
a solution to the pair of equations (|3.2p and (j3.3p . Under the condition that the equilibrium 
measure is one-cut supported, this pair co,ci yields expressions for the support and the density 
function of the equilibrium measure, but we omit the formulas. 

We note that the external field V is convex if u > 0. If n is negative, it is not but the 
construction of the equilibrium measure given above can still be carried out formally. When u is 
negative but sufficiently close to 0, we can check that the equilibrium measure constructed in this 
way is still a probability measure. When u is a large negative number, the constructed density 
function ipv{x) is negative on an interval centered at 0, and therefore not a probability density. 
This means that the external field is not one-cut regular, and our construction fails. Based 
on the analogy with matrix models without external source, the symmetry of the equilibrium 
measure and numerical simulations, we conjecture that V is one-cut regular for values of u such 
that V'y(O) > 0. 

From ()B.9p . we derive that u = — 3ci — \c\, where ci is the positive solution to (|B.9p . 
This makes u a strictly decreasing function of ci. Since cq = 0, it is easy to see that I+(0) is on 
the imaginary axis, and we denote it as I+(0) = ip {p > Qi). From the relation 

J,„eo(I+(0)) = ciI+(0) - log = 0, (B.IO) 

we derive that ci = | arctan ^ and ci is a strictly decreasing function in which means that 
It is a strictly increasing function in p. 

Like (jB.Sp . with our quartic V (using the fact that cq = 0), we have by (|3.2ip 

^(^)^|-(3cf.^ + n)log£| + 3c,s(loggj)^(loggj)'-3c?. .gC\A ^^^^^ 
I c\s^ + UCIS + 3cf s + \ s e D. 



47 



Similarly to the quadratic case, we can recover the equilibrium density using (j3.23p . In particular 
at zero we have 



^y{^) = l9M_(I+(0)) = -9 fc?I+(0)3 + (nci + 3c?)I+(0) + M 

= i (-c?/ + (1 - lb) (B.12) 

-f(l-?(P^ + i))-f(l-(^|)(a.ct.nI^)»). 

Here we used ()B.9p to pass from the first to the second line. Thus V'v(O) > if and only 
if (| + ^)(arctan j^)'^ < 1, which is equivalent to p > pc for some value pc > 0. Since u 
is an increasing function in p, this is equivalent to n > Uc, where Uc can be approximated 
numerically as Uc ~ —1.9250. Although we have not rigorously proved that for u > Uc the 
external field is one-cut regular, numerical results are convincing. When u = —1.925, the 
constructed equilibrium measure is shown in Figure [9j It suggests that around u = Uc ^ —1.925 
the transition between one-cut and two-cut equilibrium measures occurs. 




n 



Figure 9: The density function ipv of the equilibrium measure for the external field V{x) = 

4 2 2- 



C Asymptotics of pn (x) when V{x) = y 

In this appendix, we give an alternative derivation of the asymptotic results in Theorem [5] when 

2 

the external field is V{x) = The derivation is based on the contour integral formula of mul- 
tiple Hermite polynomials in Theorems 2.1 and 2.3]. This method can essentially reproduce 
all results in Theorem [2] for quadratic external field, but for brevity we only give the derivation 
for pi"^ (x) where x £ R and is away from the edges of the equilibrium measure. Although this 
contour integral method cannot be applied when the external field is not quadratic, it shows 
that the transformation Jci,co arises naturally in the uniform external source model. 

The result Theorem 2.1] states that the monic polynomial Pnix) of degree n that satisfies 

/oo 2 
Pn(x)e"'*^^e-"^dx = 0, for j = 1, 2, . . . , n, (C.l) 
-oo 

is expressed by an integral over the imaginary axis: 

P„(x) = ^ / e5(*-^')' - «*)^*- (C.2) 

y ZTTI J —too 
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When aj = as in (jl.Sp . we have, in our notations, P„ = p'n^ where V{x) = Setting 
t = s + ^ , we have 

p(r)(x) = ^ r e-^-^^-'-Us, (C.3) 

where 

x) = + i - x)2 + log(. + \- (C.4) 

2 2 2 n n 

i=l 

For s away from the interval [— ^, we have the following uniform (in s and x) asymptotic 
expansion as n — )■ oo, 



Fn{s; x) = F{s; x) + - log J ^ + O(^), (C.5) 
n ^ s - 

where 

x) = l(s + 1 - x)2 + + 1) log(. + 1) - (s - 1) log(s - ^) - 1, (C.6) 
and we take the principal branch of the logarithm and square root. Hence 

d 1 s — - 

—F{s;x) = s + - -log ^-x = 3^i{s)-x. (C.7) 

as 2 s + '2 

Below we consider the zeros s of the derivative ■^F{s;x) and express them as functions in x. 
We use the functions Ii(x),l2(a;) and their boundary values l±{x) as defined in ()1.26p - ()1.29p 
with ci = 1 and cq = ^. Note that Sq = — ^ and = as given in ()1.24p : we denote 

a = Ji,i(sa) = ^ + 2 log — - — , b = J^^i{sb) = — - — + 2 log — - — , (C.8) 

as in (j3.4p . which agree with ()B.4p . We can say the following about the zeros of ■^F{s; x): 

(1) if X > 6, then there are two zeros of ■^F{s;x): Ii(x) G (sb,oo) and 12(2;) G (|,Sf,), 

(2) if X < a, then there are two zeros of ■^F{s; x): Ii(x) G (—00, Sa) and l2{x) G {sa, —^), 

(3) if X G (a, 6), then there are two zeros of ^F(s;x): I+(x) G 71 and I_(x) G 72. 

By explicit computation, we find that for x G (—00, a) U (6, +00), then along the vertical line 
{z = Ii(x) + it \ t £ M}, attains its maximum at z = Ii(x). If we deform the contour iR 

of integration in (|C.3p to the vertical line through Ii (x) , the standard application of the saddle 
point method yields 

V27rz yii(a;)-n"^i 



_i 

n 



fixyxj+n ii / c; _|_ i 

/ ^ /i±|d,(i + o(^-i)) 

iFn(li{x);x) |-Ii(x)+n"3j / ^2 

-7^=^ / 1 exp -(s-Ii(x))2--^F(s;x) 

V27n Ji.(x)~n-ii \ 2 cfs^ 



^^nFn(li{xy,x) .Ii(x)+n 3i ^^^^ T^^^^2 
' Ii (x)— n~ 3 i 



s=Ii(a')/ V 2 



S + i 1 

f-(is(l + 0(n-2)) 



Ii(x) 



2 _ 5 



(C.9) 
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If X G (a, b), by explicit computation, we find that along the vertical line that passes through 
l-\-{x) and I-(x), ^F{z) attains its maximum at two points z = l-\-{x) and z = I-(x). (Note 
that although F{z) is discontinuous on the interval [— ^, |], ^F{z) is continuous everywhere.) 
Then we take the contour in ()C.3|1 as this vertical line. When the contour crosses the interval 
[— ^, F{z) is no longer a good approximation of Fn{z), but we can estimate the magnitude 
of Fri,{z) by other methods, and still find the vertical line suitable for saddle point analysis. The 
standard application of saddle point method yields, like (|C.9p . gives 

r- /•l±{x)+n~ii J , (rr) ^ 2. 

V" / ^n^n(s;.)rfs = e"^(i±(^);^)4^i^4^^(l + 0(n-5)), (C.IO) 



V2^i Jl±{x)-n-h Jl±(x)2 - I 



and 

\/27ri Wl+(x)-n"3j Jl_{x)-n~^i 
■,x) + 



2gn5RF(I+(x);z) 



|I+(^) 



11^ 



COS n9F(I+(x);x) + arg 




where the square roots take the principal value. It is not obvious that the asymptotic formulas 
()C.9p and (jC.lip agree with the formulas (jl.37p and (|1.4ip . To convince the reader, we show 
that ()C.9P is equivalent to ()1.37p (with k = and x G M) in the leading term. 
It is easy to check that 



Il(x)2-| 



G'o(Ii(^)) (C.12) 



where Go is the function defined in ()1.30p with ci = 1. We need also to show that F(li{x);x) = 
g(x) where g(x) is defined in (|1.33p . Since it is not hard to verify by direct computation that 
g(x) = log(x) + o(l) and F{Ii{x); x) = log(x) + o(l), we need only to show that the function 
g'(x) = G{x), defined in (jS.lOp . satisfies 

G(x) = ^F(Ii(x);x), (C.13) 
Note that by the relation x = J-,^ i (Ii(x)), we have 
F(Ii(x);x) = F(Ii(x)) 

Ii(x) + ^)log(Ii(x) + i) - (Ii(x) - ^)log(Ii(x) -I, 

(C.14) 

where we consider F 3jS Si function of n = Ii(x), and 




— F(Ii(x); x) = —F{u) — ^ = log f . (C.15) 
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On the other hand, by the identities (|3.15|) and (j3.2ip . we have 



1 



27ri Jy ^ — liix) 

di-—i , t \ di. 



By the calculation of residue, it is obvious that the first contour integral in the second line of 

()C.16p vanishes, and after some effort, we find the second contour integral has value — 27ri log ^^||^|"'"f 

Thus (|C.13p is proved, and together with (|C.12p the equivalence between (jC.Op and (jl.37p is 
obtained. 
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